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research  community  as  to  the  validity  of  these  data.  Information  on  the  correct  use  of 
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ABSTRACT 

(Distribution  Limitation  Statement  A) 

An  investigation  of  the  accuracy  of  the  Air  Force  Weapons  Laboratory 
Nonsimilar  Turbulent  Boundary  Layer  Computer  Code  through  comparisons  of  pre¬ 
dictions  with  experimental  data  is  presented.  Five  complete  data  sets  are 
selected  from  a  literature  survey  which  originally  considered  over  one  hundred 
separate  turbulent  boundary  layer  experimental  investigations.  The  five  data 
sets  include  flat  plate  supersonic  flow,  hypersonic  flow  with  and  without 
acceleration,  supersonic  flow  with  highly  nonsimilar  wall  temperature  distri¬ 
butions,  subsonic  flow  with  various  wall  blowing  rates,  and  supersonic  flow 
with  wall  blowing.  Before  the  present  study,  the  turbulent  model  in  the  AFWL 
code  was  relatively  untested  for  compressible  flows,  but  is  left  essentially 
unchanged  after  detailed  comparisons  with  other  popular  turbulent  models. 
Detailed  predictions  for  each  of  the  five  data  sets  are  carried  out,  and 
approximately  75  graphical  comparisons  are  presented  which  incluae  velocity, 
total  temperature,  and  Macn  number  profiles  at  several  axial  stations, 
plus  momentum  thickness  and  drag  coefficient  variations  along  the  flow  direction 
Overall  agreement  between  the  profile  predictions  and  the  experimental  data 
is  good  for  flow  with  and  without  blowing.  Drag  coefficient  prediction  for 
blown  flows  is  typically  below  reported  values  for  the  data  sets  selected. 
However,  these  results  are  taken  as  inconclusive  since  there  is  disagreement 
among  the  turbulent  boundary  layer  research  community  as  to  the  validity  of 
these  data.  Information  on  the  correct  use  of  the  computer  code  for  turbulent 
boundary  layer  prediction  - ,  and  on  a  new  entropy  layer  option,  is  also  presented 
The  entropy  layer  model  directly  couples  the  inviscid  flow  entropy  gradients 
caused  by  shock  curvature  with  the  boundary  layer  edge  boundary  conditions. 
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SECTION  I 
INTRODLCTION 


The  Air  Force  Weapons  Laboratory  recently  sponsored  the  development  of 
^  boundary  layer  computer  code  (reference  1)  for  use  in  predicting  multi¬ 
component,  chemically-reacting,  laminar  or  turbulent  flows  over  ablating, 
reentry  vehicle  surfaces.  This  code  has  been  designated  BLIMP  for  Boundary 
Layer  Integral  Matrix  Procedure.  The  turbulent  model  used  in  this  code  was 
recognized  to  be  relatively  untested  for  the  high  Mach  number  boundary  layers 
of  interest  to  the  Air  Force,  although  it  had  been  shown  to  be  satisfactory 
for  low  speed  unblown  and  blown  turbulent  boundary  layers  (reference  2)  .  This 
report  presents  the  results  of  a  study  to  validate  and  to  improve  that  turbulent 
model  for  compressible  and  nonsimilar  flow  cases. 

In  this  study,  the  available  turbulent  boundary  layer  literature  was 
searched  in  detail  for  experimental  data  of  sufficient  accuracy  for  detailed 
boundary  layer  profile  and  integral  parameter  comparisons.  Interviews  with 
principal  investigators  were  conducted  to  learn  of  new  data  or  unreported  data 
reduction  techniques  which  are  often  missing  from  the  usual  final  report.  Five 
data  sets  were  selected  from  approximately  one  hundred  experiments  originally 
considered. 

The  turbulent  model  in  the  BLIMP  code  was  also  re-evaluated  in  view  of 
the  numerous  papers  and  reports  on  turbulence  modeling  currently  becoming  avail¬ 
able  ir.  the  literature.  Limited  changes  in  the  existing  BLIMP  turbulent  model 
(as  described  in  reference  1)  were  made  before  the  detailed  data  comparisons 
presented  herein  were  finally  generated.  Extensive  use  of  the  BLIMP  code  on 
the  selected  test  rases  then  resulted  in  approximately  75  graphical  comparisons 
of  the  code  predictions  with  the  experimental  data.  Questions  of  code  starting 
procedure,  nodal  spacing,  etc.  were  also  addressed  in  order  to  give  a  complete 
picture  of  the  use  of  the  BLIMP  code  and  the  kind  of  accuracy  that  can  be  ex¬ 
pected  from  it. 

This  report  discusses  the  turbulent  model  study  in  essentially  the  order 
in  which  it  was  carried  out.  Section  I I  describes  the  literature  survey  phase 
including  the  criteria  used  for  selecting  the  final  five  date-,  sets.  Section 
III  describes  the  turbulent  model  used  in  BLIMP,  and  presents  some  limited  com¬ 
parisons  with  other  turbulent  models.  Section  IV  presents  the  comparisons 
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between  BLIMP  predictions  and  the  experimental  data.  Section  V  contains  some 
conclusions  and  recormendations  for  farther  work.  A  general  discussion  on  the 
use  of  the  code  for  turbulent;  fjow  problems  is  included  in  an  appendix. 
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SECTION  II 

COLLECTION  AND  EVALUATION  OF  TURBULENT  BOUNDARY  LAYER  DATA 

The  number  of  experimental  studies  of  turbulent  boundary  layer  flows 
which  have  been  conducted  since  the  turn  of  the  century  is  possibly  in  the  tens 
of  thousands ,  Of  these  numerous  studies,  possibly  1  to  10  percent  were  carried 
out  with  sufficient  accuracy  and  detail  to  yield  data  of  "acceptable"  quality. 
Of  these  high  quality  experiments,  perhaps  10  percent  include  the  type  of 
flow  configurations  and  the  necessary  instrumentation  to  provide  data  of 
interest  in  a  turbulent  model  study.  This  leaves  possibly  a  few  hundred  papers 
and  reports  which  should  be  reviewed  for  a  study  of  the  type  fceing  reported 
here.  This  survey  was  limited  (by  the  constraints  cf  both  time  and  practical¬ 
ity)  to  reports  or  papers  generally  found  or  referenced  in  the  recent  open 
literature.  Thus,  approximately  one  hundred  documents  «e.:e  surveyed  briefly 
for  pertinent  data.  Descriptions  of  test  configurations,  test  conditions, 
instrumentation,  and  the  type  and  quality  of  data  tahen  fcr  the  most  relevant, 
experiments  is  presented  later  in  this  report. 

In  the  remainder  of  this  section  of  the  report,  the  experiments  referred 
to  above  are  summarized  and  discussed  in  Section  II. 1.  Criteria  for  selecting 
the  best  data  sets  are  presented  in  Section  II. 2,  and  the  final  selected  data 
cets  are  discussed  in  some  detail  in  Section  II. 3. 

1.  SUMMARY  OF  EXPERIMENTAL  DATA  CONSIDERED 

There  are  several  ways  in  which  the  various  experiments  can  be 
described;  however,  it  is  felt  that  a  presentation  in  tabular  form  is  the  most 
efficient  and  will  be  most  useful  for  future  reference.  Thus,  two  tables  are 
presented  here  to  summarize  the  experimental  data  which  were  considered. 

Table  I  includes  flows  without  blowing,  while  Table  II  includes  flows  with 
blowing.  Within  each  table,  the  experiments  are  placed  in  alphabetical  order, 
according  to  the  principal  author's  name  in  order  to  facilitate  cross-referenc¬ 
ing  with  Tables  III  or  IV  or  the  text  of  this  report.  A  shorthand  notation  has 
been  introduced  for  some  columns  {SFB  for  skin  friction  balance,  etc.)  and  a 
"Remarks"  column  has  been  included  to  give  additional  details  about  the  experi¬ 
ment.  Often  the  remarks  refer  to  the  suitability  of  the  data  for  turbulent 
model  studies,  since  it  is  difficult  to  glean  this  information  from  the  other 
tabular  data. 
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while  the  information  provided  in  Tables  I  and  II  is  useful  for 
detailed  review,  it  doe*  not  provide  an  overall  impression  of  the  availability 
cf  certain  types  of  experimental  data  without  a  detailed  review  by  the  reader. 
Therefore,  Tables  III  and  IY  were  prepared.  Table  III  applies  to  unblown 
flows  while  Table  IV  applies  to  blown  flows.  Referring  to  Table  III,  Unblown 
“lows,  it  is  seen  that  very  few  subsonic  cases  are  listed.  This  is  entirely 
attributable  to  the  fact  that  subsonic  unblown  data  were  net  sought  in  the 
literature  search.  Matting,  et  al.  (reference  22)  plus  several  reports  of 
the  Stanford  series  present  data  on  subsonic  unblown  layers;  therefore  this 
category  was  included  for  completeness-  The  other  entry  is  from  the  1968 
ArOSP-Sttr.  ford- IFF  Turbulent  Boundary  Layer  Conference  Proceeding  (reference 
6*!.  and  is  regarded  as  top  .quality  data  useful  for  reference  purposes. 

There  is  a  noticeable  difference  between  the  number  cf  reports  on  super¬ 
sonic  flat  plate  flews  (1  <  X  <  4};  and  hypersonic  flows  (M  >4);  hypersonic 
flew*  see*  tc  have  received  xore  attention  from  the  experimentalist.  This  is 
attributable  to  the  need  over  the  last  decade  for  data  suitable  for  reentry 
calculations  afid  comparisons .  Thus,  major  test  facilities  around  the  country 
have  been  constructed  and  used  primarily  in  the  hypersonic  flight  range. 

Data  on  a — elerating  boundary  layers  were  also  gathered  during  the  liter¬ 
ature  survey  phase.  Pressure  gradient  investigations  are  often  more  unavoidable 
than  intentional  in  nary  test  facilities,  thereby  yielding  rather  uncertain  flow 
data.  However,  there  have  been  some  research  prograxs  in  which  pressure 
gradient  was  varied  systematically. 

Turbulent  boundary  layer?  with  a  ncnreactive  transparent  have  been 
investigated  fer  nary  years  as  a  possible  technique  for  heat  transfer  alleviation 
in  tucfc'ze  blades  and  reentry  vehicle  nose  tipsf  among  other  applications.  Air 
hau  been  the  most  popular  ir.jectant  for  experiments,  although,  possibly  not  the 
nest  practical  in  many  actual  cooling  systers.  The  literature  review  for  sub¬ 
sonic  blown  flows  was  r.tendec  tc  supplement  that  given  in  reference  2.  Thus 
the  Massachusetts  Institute  of  Technology  (MIT)  work  (i.e.,  references  65 
and  66)  was  not  included  for  re-evaluation.  There  is  no  shortage  of  new 
low  speed  datr  available.  The  Stanford  investigations  have  resulted  in  at 
least  13  separate  reports  or.  blown  turbulent  boundary  layers,  with  heat 
transfer  and  pressure  graciert  effects  documented  as  well.  Strangely  enough, 
the  creat  interest  in  lew  speed  bloving  investigations  has  not  carried  over 
to  the  supersonic  and  Hypersonic  flow  regions.  Particularly  at  Mach  numbers 
sreater  chan  four,  there  is  s  real  need  fer  further  basic  flow  data  suitable 
fer  turbulent  model  investigation  ar.i  development.  The  picture  is  equally 
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incomplete  in  the  nonreactive  foreign  gas  injection  category,  with  fiartle 
and  Leadon's  study  (reference  43),  where  no  profile  information  was  taken, 
being  perhaps  thf  most  exhaustive. 

Finally,  the  category  of  most  interest  for  heat  shield  ablation  or  tran¬ 
spiration  applications  is  the  reactive  blowing  data.  Work  is  just  getting  under¬ 
way  in  low  speed  flows  with  the  University  of  Utah  duplication  (references  41, 

49,  and  59)  of  the  original  Wooldrige  and  Muzzy  experiments  (reference  6:). 

There  is  little  basic  data  of  interest  in  supersonic  or  hypersonic  flows. 

2.  CRITERIA  USED  FOR  DATA  SELECTION 

The  criteria  used  for  judging  the  usefulness  of  experimental  data  must 
of  course  depend  on  what  use  is  to  be  made  of  the  data.  For  turbulent  model 
studies,  it  is  of  interest  to  model  the  details  of  the  flow  through  the  entire 
boundary  laver,  therefore  profiles  of  the  important  boundary  layer  variables 
become  of  primary  importance.  In  the  usual  experiment,  state-of-the-art  in¬ 
strumentation  is  used  to  measure  stagnation  pressure  and  stagnation  temperature 
behind  a  normal  shock  (created  by  the  probe  itself)  for  supersonic  flows, 
although  for  adiabatic  wall  conditions  local  total  temperature  can  be  assumed 
equal  to  the  reservoir  temperature  With  only  a  small  error.  These  pressure 
and  temperature  data  are  used  to  arrive  at  other  more  fundamental  variables 
such  as  velocity  and  temperature  or  enthalpy.  Thus,  data  selected  for  these 
turbulent  model  studies  must  at  least  include  measured  profiles  of  one  or  more 
of  the  important  boundary  layer  variables 

Other  criteria  on  which  to  judge  the  data  are  more  subtle  anc  can  only 
be  described  in  a  fairly  general  way.  Herein  is  presented  a  list  of  items 
that  have  been  considered  in  viewing  a  set  of  experimental  data. 

•  The  wind  tunnel  should  provide  a  shock-free,  low  freestream 
turbulence  isentopic  flow. 

•  The  boundary  layer  flow  upstream  of  the  test  section  shou^i 

be  well  characterized  if  it  eventually  forms  the  boundary  layer 
for  testing  (as  in  wind  tunnel  wall  boundary  layer  testing)  . 

•  The  test  section  or  model  should  be  of  known  flatness  and  1 
hydraulically  smooth. 

•  Leading  or  trailing  edge  effects  should  be  examined  to  determine 
their  influence. 

•  The  two-dimensionality  of  the  model  -  wind  tunnel  -  test  section 
combination  should  be  considered. 
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•  The  flow  should  be  fully  turbulent  at  the  test  section,  free  of 
any  transitional  or  boundary  layer  trip  effects. 

•  Where  transpiration  effects  are  investigated,  the  injection  must  be 
uniform  or  known  to  within  a  few  percent. 

•  The  displacement  effects  of  the  boundary  layer  flow  on  the  inviscid 
flow  must  be  completely  known  -  particularly  with  injection. 

•  Probe  sizes  must  be  small  compared  to  boundary  layer  sizes. 

•  Disturbance  of  the  flow  at  the  wall  by  wall  instrumentation  must 
be  held  to  a  minimum. 

•  The  data  should  appear  relatively  smooth  and  consistent  with  a 
minimum  of  scatter. 

•  Where  possible,  the  integral  conservation  equations  should  be 
checked  and  satisfied  experimentally  to  a  reasonable  percent 
error. 

The  points  mentioned  above  are  essentially  the  desired  features  of  an 
acceptable  boundary  layer  experiment.  No  attempt  was  made,  however,  to 
develop  a  rating  system  for  the  various  data  sets.  Rather,  an  overall 
impression  gathered  from  reading  about  each  experiment  was  formed,  notes 
were  made  regarding  each  data  set  where  necessary,  and  in  many  cases, 
telephone  or  personal  interviews  were  held  with  principal  investigators. 
Important  features  of  each  data  set  have  been  noted  in  the  "Remarks"  column 
of  Tables  I  and  II. 

3.  SELECTED  DATA  SETS 

A  minimum  of  five  data  sets  was  desired  for  turbulent  model  comparisons . 
The  approach  used  in  selecting  these  data  sets  was  an  attempt  to  provide  data 
which  test  the  turbulent  model  in  a  variety  of  flow  situations  which  are 
encountered  in  flight.  This  obviously  means  that  comparison*;  could  not  be 
carried  out  for  each  category  of  Tables  III  and  IV  since  there  are  mere  than 
five  categories  containing  good  data.  The  five  data  sets  selected  were: 

1.  Coles'  JPL  experiments  (reference  9  -  Supersonic  flow  with  zero 
pressure  gradient. 

2.  NOL  experiments  by  Lee.  et  al.  (reference  17)  and  Brott,  et  al. 
(reference  6)  -  Hypersonic  flow  with  and  without  pressure  gradient. 
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3.  TRW  expe  i Bents  run  at  AEDC  (reference  18. '  -  Hypersonic  flow  with 
highly  nonsimilar  wall  conditions. 

4.  Stanford  experiments  by  Simpson,  et  al.  (reference  331  -  Low  speed 
flow  with  nonreactive  blowing,  including  steps  in  blowing. 

5.  Cambridge  experiments  by  JeroauLn  (reference  48)  -  Supersonic  flow 
with  nonrcactive  blowing. 

Detail;  tn  the  five  selected  data  sets  are  presented  in  the  reaaind  ;r  of  this 
section.  The  descriptions  below  are  supplemented  by  three  more  tabular  compari¬ 
sons.  Table  V  presents  the  particular  test  conditions  for  each  test  run  used 
for  comparisons.  Table  VI  presents  further  information  on  the  instrumentation 
used  in  each  of  the  five  data  sets,  while  Table  VII  compares  the  data  reduction 
techniques . 

a.  Coles  -  Supersonic  Flow  with  Zero  Pressure  Gradient 

Coles'  data  (reference  8)  at  two  Mach  numbers,  M  -  2.6  and  M  «  3.?, 
were  selected  as  baseline  cases.  Coles'  experiments  were  carried  out  in  the 
Jet  Propulsion  Laboratory  20-inch  supersonic  wind  tunnel,  whose  flexible  nozzle 
walls  were  adjusted  to  give  a  nearly  constant  pressure  flow  over  the  adiabatic 
flat  plate  model.  The  test  model  consisted  of  a  30-inch  long  flat  plate  sharp¬ 
ened  with  a  half -wedge  on  the  upper  surface  at  each  end;  the  lower  test  surface 
was  flat  over  the  entire  length.  This  plate  was  inserted  into  the  tunuel,  lo¬ 
cated  on  the  centerline  and  pitched  0.17  degrees  nose  down.  Various  types  of 
boundary  layer  trips  were  employed  by  Coles.  For  the  tvo  comparison  runs  cited 
herein,  the  trip,  or  fence,  consisted  of  0.014-inch  diameter  vires  spaced  1/4 
inch  apart  which  projected  about  0.1  inch  beyond  the  leading  edge  of  the  plate. 

This  experiment  emphasized  accurate  evaluation  of  the  skin  friction  coef¬ 
ficient  with  a  floating  element  gage.  Only  total  pressure  profiles  were  meas¬ 
ured:  total  temperature  was  assumed  to  be  constant  to  determine  velocity  pro¬ 
files  across  the  boundary  layer  for  the  adiabatic  plate,  as  described  in  Tables 
VI  and  VII.  Preestream  conditions  were  assumed  based  on  isentropic  expansion 
of  a  y  =  1.4  fluid.  Profile  measurements  were  taken  at  only  one  station  21.5 
inches  from  the  leading  edge.  Skin  friction  gages  were  placed  at  5.5,  13 
and  24  inches.  Static  pressure  taps  in  the  plate  were  used  to  determine  axial 
pressure  variations. 
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b .  Naval  Ordnance  Laboratory  -  Hypersonic  Flow  Without  and  With 

Pressure  Gradient 

Many  reports  and  documents  were  found  which  presented  hypersonic  flow 
data,  and  most  of  them  included  profile  measurements.  After  discussions  with 
various  investigators  in  this  field,  the  recent  zero  pressure  gradient  experi¬ 
ments  at  the  Naval  Ordnance  Laboratory  by  Lee  et  al.  (reference  17  j  and  the  ex¬ 
tension  to  pressure  gradient  flows  by  Brott  et  al.  (reference  6)  were  selected. 
Both  of  these  investigations  were  carried  out  using  the  NOL  Mach  5  boundary 
layer  channel,  which  is  a  two-dimensional  half  nozzle  with  a  flexible  wall 
forming  the  nozzle.  Boundary  layer  measurements  are  carried  out  on  the  opposite, 
flat  wall ,  with  moderate  heat  transfer  to  the  wall  controlled  by  circulating 
wall  coolant.  This  stainless  steel  test  surface  is  eight  feet  long  and  tapered 
from  12  to  13.5  inches  wide.  Lee’s  tests  were  carried  out  at  four  stations  be¬ 
tween  48  and  92  inches  downstream  of  the  throat,  giving  naturally  turbulent 
boundary  layers  from  two  to  four  inches  thick.  Brott's  measurements  were  at 
five  stations  between  50  and  84  inches,  giving  boundary  layer  thicknesses  of  1.3 
to  3.9  inches. 

A  wealth  of  instrumentation  was  used  in  these  NOL  investigations.  Profile 
•  instruments  included  a  total  pressure  probe,  a  static  pressure  probe,  an  equi¬ 
librium  temperature  probe  (measures  local  recovery  temperature),  and,  in  Brott's 
investigation,  a  fine  wire  thermocouple  probe.  Wall  instrumentation  included 
water-cooled  skin  friction  balances  upstream  of  each  profile  measurement  location, 
and  wall  heat  flux  gages. 

c.  TRW  -  Hypersonic  Flow  with  Step  in  Nail  Tareperature 

The  Arnold  Engineering  Development  Center  (AEDC)  Supersonic  Kind 
Tunnel  “A"  with  a  4C-inch  by  40-inch  test  section  provided  the  main  air  flow 
for  the  TRW  experiments.  This  tunnel  was  operated  in  a  continuous,  closed-circuit 
mode.  The  nozzle  used  was  atu  automatically  driven,  f lexible-plate-type. 

Some  cf  the  test  details  are  contained  in  an  AEDC  report,  reference  67. 

The  test  node!  was  a  20.09-ir.ch  I  D.,  49 . 48- inch  long  hollow  cylinder 
aligned  with  the  f low  to  +  0.05  degrees  _n  pitch  and  ya>  about  the  tunnel  center- 
line  aod  supported  by  two  struts  attached  to  the  downstream  end.  In  cross  sec¬ 
tion,  the  leading  edge  was  a  half  wedge  on  the  outer  surface.  The  inner  (test) 
surface  was  flat  along  its  entire  length  with  a  surface  finish  of  100  --inches. 
Seven  separately  controlled  cooling  (heating)  copper  coil  circuits  were  counter- 
wrapped  or.  the  outside  surface.  The  409°F  step  in  wall  temperature  was  made  be¬ 
tween  23.0  and  24.0  inches  tros  the  leading  edge  with  four  of  the  seven  circuits 
forward  of  the  step.  Liquid  nitrogen  was  circulated  ir.  the  cold  wall  region; 
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pre-conditioned  water  was  used  to  achieve  near-adiabatic  conditions  in  the  *fcot” 
region.  For  the  temperature  step  tests,  the  circuits  in  each  region  were  naci- 
folded  together.  A  0.375-inch  wide  by  0.006-incJ  thick  piece  of  fiberglass 
tape,  serrated  at  its  leading  edge  and  located  on  *■»*“  inner  surface  approximately 
0.125  inches  from  the  cylinder  leading  edge  provided  a  boundary  layer  trip. 

Profile  instrumentation  for  these  tests  included  a  total  pressure  probe, 
static  pressure  probe,  shielded  stagnation  temperature  probe,  and  a  constant 
current  hot  wire  anemometer.  The  total  temperature  probe  size  '0.06C  inch  in 
diameter)  precluded  measurements  very  near  the  wall.  Wall  condition*  were 
measured  with  static  pressure  transducers,  surface  temperature  thermocouples, 
and  Garden  heat  transfer  gages  for  wail  heat  flur.  Skin  friction  was  measured 
with  38  individual  Stanton  tubes. 

d.  Stanford  -  Zero  Pressure  Gradient,  Subsonic  Flow  With  and 
without  Uniform  Slowing,  ar.d  Witn  a  Step  in  Blowing 

The  heat  and  mass  transfer  apparatus  at  Stanford  university ‘s 
Mechanical  Kngineenng  Department,  a  facility  built  specifically  for  high  acur- 
racy  subsonic  boundary  layer  experiments  with  blowing  and  sucticus.  was  used  by 
Simpson  (reference  33?  for  the  selected  uniform  blowing  tests.  A  complete  des¬ 
cription  of  this  apparatus  is  given  bv  Moffat  (reference  25? .  The  nain  flew 
system  was  a  continuous  operation  open-circuit  consisting  of  an  inlet  filter, 
flow  control  valve,  blower,  heat  exchanger,  screens  (for  flow  straightenir • . . 
settling  plenum  chamber,  primary  nozzle,  and  test  section.  The  rectangular  test 
section  was  foraed  of  two  fixed  side  walls,  a  fixed  flat  bottom  wall  which 
included  the  porous  test  surface  and  an  adjustable  top  wall  cf  plexiclas  con¬ 
trolled  by  a  bank  of  jacks. 

The  test  surface  consisted  of  24  individual,  I '4 -inch  thick,  porous  place 
sections  which  were  10  inches  wide  by  4  inches  long,  resulting  ;o  a  total  length 
of  95  inches.  The  plates  were  constructed  of  C . 005- inch  sintered  bronze  spheres 
which  provided  a  surface  finish  of  20C  u-inches.  Final  assembly  cf  the  plates 
left  a  0.032-inch  gap  between  adjacent  sections.  This  gap  was  filled  and 
smoothed  br.t  was  not  porous. 

The  secondary,  or  injection,  flow  system  consisted  cf  an  inlet  filter, 
blower,  heat  exchanger,  and  f lew  header  followed  by  24  ~ serrate  flow  control 
valves,  flow  meters  and  delivery  tubes,  one  set  for  each  plate.  A  porous  sin¬ 
tered  bronze  pre-plate  was  used  tc  obtain  mifoia  pressure  in  the  small  plenum 
behind  each  porous  surface  plate.  Localized  transpired  flow  rates  were  found 
to  be  *  c  percent  ever  the  center  6  inches  cf  the  spar..  The  heat  exchangers  in 
the  rain  and  secondary  systems  were  controlled  to  provide  equal  flow  temperatures 
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A  trip  of  ]/>-■ .act  wide  eocic  srid  cadorsdaa  pacer  Ix«e«<f 

at  f\*  exit  tdx  cf  the  tnetitix  s*cti»  20:21s?  tie  rcrrle  tc  a*  tes : 
section.  9k  tniiia;  eiy  of  tie  J.IU-ikS  title*  trip  «»  L  t-i* Ji  s?sv»bi 
ci  a*  gawa  test  Mctix. 

Omidenblc  tap&asis  was  placed  «  saig ns?  expcnxeatal  errors  is 
itt  castrciiaa  of  tie  Suafori  particxlarly  is  grc-ridisr  s.fxa 

isjertioa  a ass  flew,  icowttiag  for  all  scarce?  or  losses  of  tcarcy.  a»£  is 
prowidu?  «  lew  fcerbwlenoe  anfen  sai s  flaw.  flow  rates  and  slate  tejper atares 
were  isdindullr  searered  for  eaca  of  tie  21  pc.rtw  swoests.  Frews  ere  profiles 
were  tabes  «tt  a  flattened.  -di*-ireft.  by  .i3>io±  probe  at  fane  axial  loca¬ 
tions  for  eaci  test,  dsris?  »t: ct  tie  rreestresni  ccai;t:cci  were  wustaired 
eirtaally  ocAstaat.  Static  taptratem  were  aeasrred  wits,  at  exposed  tsersc- 
coaple  probe  aliased  with  tse  flow.  Static  pressure  taps  located  alecs?  tie 
side  walls  were  used  ta  obta.tr  tie  desired  asifcrs  flaw  err  tree  of  tie  ascer  scr- 
face.  2ee  to  tie  cse  of  pcrccs  plates,  sils  friction.  was  ret  aeaxcrsd  directly. 

was  copatad  freat  sae  ~nr>w » -.m  irteprxl  equation  aud  by  a  tiskot  sublayer 
setioc  fsee  refereoce  53'  . 


e -  derate  -  Sape rsoric  claw  Kiri  aad  eiticst  Sscfccs  Jlcwrg? 

-Jemair's  coper iaw-tts  were  carried  cat  is  tie  la abridge  S2irersity 
Sxriaecrla?  Depsrtawat  Ssoarsacnc  Xiad  Tenrel.  as  irterxi  treat  iliwdcwo  tsr&sl 
w-iti  dost  W  seconds  of  useful  t*s*  tise.  Tie  aoxc le  was  twcsSiaece i anal  with 
fixed  side  aad  tsx»  serf  aces.  Tie  tap  surface  was  contoured  bat  act  adjustable; 
tl<  Icwr  surface  wfeicA  costajsed  tie  first  apgt*d  test  surface  was  flit  sac 
could  be  adjusted  vertical  It  is  a  pitxii  xde.  Tie  rectargrler  test  section 
was  *.5  iacbcs  be  4.5  isefee' 


ferasus  used  two  tea  .-faces;  a  solid  plate  f cr  tie  zero  riewtr?  casts 
*  l.'4-isrft  it*  ex  relied,  *cdrosist « *  sir  ter  ed  erma*  plate  wed*  iron  2.5 
iu.croe.  softer es  foe  blower?  tests.  Ti.~ese  plates  were  ;.3J  ir«S*s  wide  iy  12. € 
ireftes  ~.coc  aad  were  inserted  irtc  tie  tettem  lire  rf  tee  sczxle  wits  tie  Lead- 
13?  edge  3.4  incises  dewnstrexat  cf  tie  nxcle  tiroac-  Tie  jirctlsc.  ietweer-  tie 
aotcie  wall  ard  test  serf  ace  was  fi.led  wits  paste  *?<•  ccctad  flat.  A  oc^sadary 
layer  vrip  was  sot  eeqslcyei. 


Injection  air  was  tabes  free  tie  •***  -  tmrel  sorely  passed  tirosari  a 
flow-  control  t*1t*  m»*  a  xeterir?  rnfo.ee.  Two  sorts  carried  tie  flow  ir.ee 
separate  pie aat  cfSarrbiers  erstreas?  or  a  octets  pre— place.  A  sircle  sierra  is 
ere  itber  side  was  ccswk  tc  tie  ertire  test  serf  ace  plate.  Kass  floss  rariai.ocss 
area?  2.5-onci  cicxeter  sauries  of  tie  test  sxrfaoe  raterial  ret  tie  test  ser¬ 
iate  itself’  were  *  5  yerrert .  to  local  variant  reas  ireaerts  of  tee  trtcal 
test  slate  were  :*ccrtcc. 
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a  p:tet  probe  as d  an  eqiiUtrias  cone  recovery  temperature 
probe  'similar  to  tbe  SX  dear “i  to  ofctai-  tre  booadtry  layer  profiles.  After 
several  wwareaeats  Aids  sbcwcg  eacellec:  a^reeaeot  vith  tire  static  temperature 
relation 


total  temperature  profiles  «ro  sot  wasurci.  Ve'ooity  and  temperature  profiles 
*«r<  ietenuaed  iteratively  os: oc  the  above  relation  and  the  Xach  cudber  profiles. 
rt:t  fr.cuce  was  act  seescreo  directly,  bet  calculated  fro*  the  aranentuar  inte¬ 
gral  e^atict.  the  Sooser  aod  Short  T’  netbed.  and  a  transformation  seethed  des¬ 
cribed  by  Jeroaia  "reference  If  ' .  The  rass  flcv  through  the  entire  plate  ves 
Teas need  vith  a  tisrlt  calibrated  ctifice. 
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SECTION  III 


TURBULENT  MODEL  SELECTION 

The  BLIMP  code  solves  the  partial  differential  equations  expressing  the 
conservation  of  mass,  momentum,  energy,  and  species  for  a  viscous  fluid.  In 
turbulent  flows,  the  equations  are  written  with  velocity,  enthalpy,  end  species 
concentrations  expressed  in  terras  of  mean  and  fluctuating  components,  i.e., 

u  =  u  +  u’  (2) 

The  time  averaged  equations  of  motion  then  include  time  averaged  products  of 
fluctuating  quantities  such  as  u’2#  which  distinguish  the  equations  of  motion 
from  their  laminar  flow  counterparts.  Turbulent  modeling  involves  the  mathe¬ 
matical  description  of  these  new  flow  variables  such  that  the  number  of  govern¬ 
ing  equations  Hatches  the  number  of  unknowns. 

The  BLIMP  code  has  been  written  to  accept  eddy  viscosity  or  mixing  length 
descriptions  of  the  turbulent  fluctuation  terms.  The  eddy  viscosity  approach 
characterizes  the  turbulent  velocity  fluctuation  as  being  related  to  local 
velocity  gradient: 

-  hv)  ‘u‘  £  pen  |~r  {3) 

Mixing  length  is  related  to  eddy  viscosity  by  the  relation 

e  =  |H  (4) 

ra  ?y 


Many  eddy  viscosity  and./ or  mixing  length  models  have  been  reported  in  the  recent 
literature  (references  68  to  74)  .  in  this  section  of  the  report  three  widely 
used  caodels  axe  reviewed  briefly  and  a  "best"  combination  ot  eddy  viscosity 
and  mixing  length  relations  and  associated  constants  is  chosen  for  the  data 
ccBoarisons  of  Section  IV. 
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1.  DESCRIPTION  OF  LEADING  TURBl’LENT  MODELS 

a.  Aerotherm  Model  (reference  1) 

The  Aerotherm  turbulent  model  was  formulated  specifically  for  bound¬ 
ary  layer  flows  with  blowing.  It  was  first  reported  in  reference  2,  and  was 
limited  to  incompressible  flows.  Turbulence  is  described  in  the  wall  region 
by  a  mixing  length  equation,  which  approaches  Prandtl's  expression,  I  =  k^ , 
far  from  the  wall.  The  incompressible  expression  is 


£* 

Ja 


(5) 


The  constant  k  is  taken  as  0.44,  while  y  ,  related  to  a  laminar  sublayer  thick- 
m  a 

ness,  is  held  constant  at  11.823.  Shear  stress  x  is  taken  to  be' the  local  value 

rather  than  the  wall  stress,  x  .  This  expression  allows  a  smooth  transition 

from  zero  turbulence  at  the  wall  to  large  turbulence  in  the  "law  of  the  wall" 

region  of  the  boundary  layer.  The  mixing  length  distribution  is  automatically 

changed  for  blown  flows  due  to  the  blowing  effect  on  shear  stress.  Constants 

v+  and  k  remain  the  same  with  or  without  blowing. 

Ja  m 

For  compressible  flows,  equation  (5)  was  changed  to  deal  with 
the  products  pi  and  oy: 


{€) 


Constants  y*  and  km  are  left  at  their  incompressible  values  cf  11.823  and  0.44, 
respectively.  As  with  incompressible  flows,  mixing  length  distribution  is 
changed  for  blown  flows  implicitly,  by  the  influence  of  clewing  on  looal  shear 
stress,  T. 

The  equation  for  mixing  length  is  solved  and  then  used  to  find  eddy 
viscosity  from  equation  (4).  This  wall  region  eddy  viscosity  is  used  fro r. 
the  wall  on  out  into  the  boundary  layer  until  it  exceeds  the  wake  region  eddy 
viscosity  expression 


t. 

E 


0.018  u  i* 
e  x 


\7) 
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at  which  point  the  constant  wake  value  is  adopted  fcr  the  remainder  cf  me 
boundary  layer.  In  equation  {7*.  is  the  local  boundary  layer  edge 
velocity  and  5?  is  the  kinenatic  cisplacener.t  thickness. 


* 


In  the  energy  and  species  equations,  terns  involving  tine  averaged 
fluctuations  of  enthalpy  and  species  concentrations  appear  and  are  related  to 
the  mean  flew  variables  by  ccliniticns  of  eddy  conductivity  and  eddv 
diffusivity: 


=  - 


In  the  Aerotherx  sod*!,  eddy  ccndectivity  and  ciiffusiv:  ty  are  siicclr  r-~I>ted 
tc  eddy  viscosity  by  constant  turbulent  Prarsdt!  arc  Sitridt  rrfetis. 


A.  value  of  C.75  fer  each  cf  these  vas  reccroe.-~.iec  in  reference  I. 

The  .lerotkeir  rod  *  1  has  ret  had  tne  men  literature  eiroscre  that  the 
ether  two  models  considered  here  have  had.  Tnercfcre.  in  erder  to  verifv  it= 
usefulness,  soma  results  fr-er  previously  run  precrarr  r-redi mens  are  .rclvded 
here.  Figures  I  to  1  creser.t  cceciriscnv  with  flat  plate,  uniicvr, 
incompressible  flew  data  taker,  by  Jfiechirct  and  Tillran  reference  is;  . 
fieghardt  and  ?ilbtar,  case  was  inc.uded  in  l?6s  Ar^S?.-?tar f cr d- 
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ouity  Ratio  Profile* 

ir/hardt,#  Cubnonlu  flow,  aero  Prueaure  Gradient 


Linear  Los  Velocity  Ratio  Profile 

Wieghardt ,  Subsonic  Plot/,  Ze  :o  rressuro  Gradient 


ri<ni'rf  1.  TM«k «*>«»•  KmvhoI'J  '  v«  ,  »'*•< *• 

Wi  ti'jiifli'ilt )  inh#i»iili»  Klfiw,  r«r«  I'ramwr* 
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Turbulent  Boundary  Layer  Conference/  and  is  regarded  as  a  top  quality  flat 
plate  experiment.  The  excellent  agreement  between  the  Aerotherm  model  and 
these  data,  both  In  profile  and  drag  prediction/  is  apparent  from  these  figures. 
Comparisons  with  subsonic,  blown  boundary  layer  data  taken  at  MIT  by  Kendall 
(reference  66}  are  shown  in  Figures  5  to  8.  These  data  were  essentially 
the  basis  for  the  Aerotherm  wall  law  model,  although  comparisons  with  other 
data  have  been  shown  to  be  equally  good.  As  with  the  unblown  case,  profile 
comparisons  are  excellent.  Drag  is  underpredicted  by  about  10-20%  in  com¬ 
parison  to  the  reported  C^/2  values  for  this  case;  however  the  prediction  is 
well  within  the  error  bands  associated  with  these  values. 


b.  The  Cebeci  -  Smith  Model  (reference  68) 

In  the  Cebeci-Smith  model,  the  boundary  layer  is  also  characterized 
by  inner  and  outer  regions.  In  the  inner  region,  a  mixing  length  approach  is 
used,  based  on  the  Prandtl  mixing  length  law  as  modified  by  Van  Driest 
(reference  75) .  Cebeci  and  coworkers  have  extended  and  modified  the  basic 
Van  Driest  law  to  include  pressure  gradient,  blowing,  and  variable  turbulent 
Prandtl  number  effects.  The  inner  region  mixing  length  expression  is 


i  •  kmy[l  "  «xp(-y/A)] 


(13) 


(14) 


(15) 


(16) 


30 


31 


Figure  5.  Velocity  Ratio  Profiles 

Xendall,  Subsonic  Flow,  F  »  0.00325 
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Figure  7.  Momentum  Thickness  Reynolds  Number  vs.  Streamwise  Location 
Kendall,  Subsonic  Flow,  F  =  0.00325 


-  Ftadtl 


.  UtfS  M  1.4 


at 


is  Ltu 


is  nflacad  is  tte  «ti«  off  s.  is  ur  efftcti  tie 

Cscbsss if  tte  sntt  cf  Untie;  dtcmsts  t&ids  reselti 
sacs  npiif  igw*e»  ts  t&t  Msttl  si  sis;  sjfits- 


Utkf  s 

Cfottiflt  I?)  IS 


rsjxas.  s 


tddr  siscos icy 


sisilsr  Cc 


Si  .  t.SIM 


T&is  «Mr  viscosity  is  «n?tipli>i  fey 


*s  Jrsfsrescs  Til  is*  tm  tree cy 


y-fi-s.s(*V] 

*  V  *  /  t 


ccyscst  * 
tfciefe  seise. 


odor?  lift  ttiritom.  as  is  t±*  Arrodfesxs  sod*!.-  ti>e  i 
os  is  oood  satil  it  »irs  ‘i  cSs  cstsx  rsyioo  *j» 


•ea.  *c 


to  Vis 


coBCactiTitf  is  ’am?  stor  as  exsstsscisa.  o£ 
r§f  »fiU«t  nbfk.c  to  as.  oodlistscy  bocadttr 
las  oi  eddy  viscosity.  Tfefs  asslysis 
sst  Pxsodtl  lister 


**- 


-  MpHAil 

-  «jqp£-wssy3} ; 


=r  jh  U/u..  -~V  -~V 

i !  n  v  vs)  j  l  m 


^  «  is  energy  siring  length  csatist ,  tun  u  8. *4 


?r  a  sslacolir  Frudtl 


35 


Six*  lebeei  docs  sat  coesidcr  -ornSr—rKycncora  flows,  be  docs  cot  ^rt  a  sodcl 

for  rcrbtticssr  feirif:  imsttr. 

Tfe*  Ctfctci-iaitv  tartalwt  node!  nliditr  is  well  doc  —  nil  u  softer- 
cos  opes.  I  iters  tore  publications  Preferences  77  to  79). 

c.  9esha*il  sad  Itebdif  Model  {reference  19)  . 

Ms5a*ll  is d  Sretwiti  esc  a  Kxiay  length  rode  I  foe*  ft*  «t:n  hxafirr 
lacer.  Se*r  tie  well  ;to  y.-  j  *  S-15.  Pnafti’t  aiii??  leayti  torts* :cc 
aotifie?  by  Tea  Driest'*  dspia;  facties  is  used: 

i  *  *.8yfl  -  esp'-y/Xb}  up: 


A 


local  sad  density  ere  iacladef  retier  rSa-  wall  ralaes.  The  X*  Talc*  is 

ttu*  fro*  a  fit  cf  nyerieuul  data  ae  blows  -  tdscsic  bocadary  1 mu,  ts 
socmb  ts  t-  tic  eijetiaeiul  date  are  those  of  Smpsoc..  et  al.  iref- 

eresoe  3 1’*  sad  tie  data  reported  rr  refertcc*  2. 


la  th*  far  wall  reject  iy/f  »  i.J’ .  sixiaj  length  is  iii—f  to  be  i 
fraction  c:  acosdary  layer  thic^ress  sad  the  iacosprersible  focst  tartor,  =*  * 
**« 


^  v> 


figure  1?  slews  the  curse  used  to  dstentire  i/i-  this  cmre  is  tie  "sraadratic 
vsriatica"  of  reference  69,  vbici  has  beer  adcctec  as  tie  standard  aooel 
at  XstSA/Iatsaley  «'  reference  85?.  rt*  5f  ^  I/i  variation  as  a  forctioc  of 
iacospressibie  for*,  factor  is  as  ittersst  to  account  for  sooeq-ailibrius!*  rf  fecit 
i-  she  profile  shape.  Thus,  it  is  o-nst  useful  in  strongly  sossiailar  flow  situ- 
iticas. 

A  straight  I  iae  seasert  for  i/i  is  used  to  join  the  inner  and  far  wall 
recicas  between  y/i  *  0-1  a ad  y/3  =  C.3. 

Steeru:lisr:sjs  as  used  here  refers  tc  the  fluid  aachanical  state  rather 
than  chenical  state . 
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Figure  10.  Variation  of  Mixing  Length  for 
Bushnell-oecVwith  Wake  Model 
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Bothntll  «U  Beckwith  use  a  constant  turbulent  Prandtl  number  but  case 
it  on  total  enthalpy  fluctuations  rather  than  static  enthalpy  fluctuations. 
That  is. 


Since  no  species  equations  are  solved,  no  turbulent  Scheldt  number  is  specified. 

the  validity  of  the  Busbnell-Beckwith  model  under  several  different 
kinds  of  flow  condi tior*  hr.  been  presented  in  reference  69. 

2.  TURBULENT  JtODEL  COMPARISONS  FOR  THE  LAW  OF  THE  WALL  REGION 

a.  NALAW  Program 


In  order  to  examine  the  differences  between  the  near  wall  mixing 
length  formulations ,  a  small  computer  code  (WALAW)  was  Written  which  solves  the 
one-dimensional  continuity,  momentum,  and  energy  equations.  These  equations, 
written  for  the  wall  region  where  streomvise  changes  jure  small  compared  to 
changes  normal  to  the  wall,  are 

continuity 


pv  =  pwvw 


momentum 


Pwvwu  -  p(u  +  Cn) 


nr  3y 


A  perfect,  ideal  gas  equation  of  state  is  assumed,  and  Runge-Kutta  integration 
is  used.  Inputs  to  the  code  are  tw,  qw,  Tw,  and  (pv)w.  Also,  a  law  for  the 
eddy  viscosity  variation  must  be  included  in  the  set  of  equations.  The  coda 
includes  tne  three  wall  region  eddy  viscosity  relations  described  above,  howt-ver 
a  constant  turbulent  Prandtl  number  equal  to  0.9  has  been  assumed  in  all  three 
models. 
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b.  Results  for  Three  Near  Wall  Models 

The  differences  between  the  three  near-wall  mixing  length  models  are 
exemplified  by  the  curves  of  Figures  11  through  13.  The  WALAW  code  was  run  for 
a  50  ft/sec,  adiabatic  flat  plate  flow  with  three  different  blowing  rates, 
i3w'\i/peue  31  0.004,  and  0.010.  The  turbulent  models  in  these  figures  are 

labeled  as  follows: 

A  =  Aerotherm  model 
B  =  Bushnell -Beckwith  model 
C  =*  Oebeci-Smith  model 

In  Figure  11  where  no  blowing  is  included  it  is  seen  that  all  three  models  result 
in  nearly  the  same  mixing  length  variations.  The  Prandtl  mixing  length  equation 
is  essentially  matched  at  about  y+  *  100.  In  Figure  12  where  a  moderate-to- 
ctre>'g  blowing  rate  of  0.004  has  been  input,  all  three  i+  curves  h^ve  shifted 
toward  the  Prandtl  law,  but  the  Bushnell  and  Cebeci  models  show  more  shift  than 
the  Aerotherm  model.  In  Figure  13  where  a  very  strong  blowing  rate  of  0.010 
is  included,  Models  B  and  C  have  adopted  the  Prandtl  line,  1  =  kmy,  over  the 
entire  range  of  interest.  Model  A  shows  .  similar  trend  with  blowing,  but  does 
not  respond  as  jvadily.  Thus,  it  can  be  expected  from  these  results  that,  under 
no  blowing  conditions,  there  will  be  very  little  difference  between  near  wall 
predictions  with  the  three  milling  length  models.  As  blowing  increases,  differ¬ 
ences  will  become  apparent.  When  constrained  to  the  same  wail  shear  value, 
the  Aerotherm  law  will  give  lower  turbulence  levels  which  in  turn  result  in 
largar  gradients  (3u/3y,  3T/3y  etc.)  than  the  Cebeci  and  Bushnell  models. 

These  conclusions  about  the  similarities  between  predictions  for  unblown 
flows  are  borne  out  by  the  WALAW  results  of  Figure  14  through  19.  Figure  14 
includes  a  velocity  profile  comparison  with  the  Kieqhardt  and  Tillman  data 
(reference  38).  All  three  models  show  good  agreement  in  the  "law  of  the  wall" 
region,  which  extends  to  about  0.5  inch  for  the  station  selected.  Beyond 
that  point,  the  one-dimensional  equations  of  motion  are  no  longer  valid. 

Figure  15  includes  comparisons  with  Coles  (reference  8)  Mach  3.7  adiabatic 
flat  plate  data,  which  shows  how  well  these  laws  apply  to  compressible  flows. 

At  lvast  for  adiabatic  flows,  any  of  the  three  models  appears  to  be  satis¬ 
factory  in  the  supersonic  range. 

For  hypersoni  •  flows,  the  experimental  data  uncertainties  become  greater 
since  the  flow  _  self  is  more  difficult  to  work  with.  Figure  16  includes  predic¬ 
tions  from  the  three  theoretical  models  compared  to  an  experimental  Mach  number 
profile  measured  by  Lee,  et  al.  (reference  17).  Mach  number  comparisons  are 
shown  since  this  information  can  be  calculated  directly  from  the  pitot  tube 
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Figure  11 .  Mixing  Length  Variation 


Figure  12.  Mixing  Length  Variation  With  F  =  0.004 
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Figure  13.  Mixing  Length  Variation  With  F  =  0.010 
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Figure  16.  Linear-Log  Mach  Number  Profile 

HOL,  Lee,  et  al..  Run  12902,  Hypersonic,  Zero  Pressure  Gradient  Flow 
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Figure  18.  Linear -Log  Mach  Number  Frofile 

NC'L,  Lee  et  al . ,  Run  6211,  Hypersonic,  Zero  ?ressure  Gradient  Flow 


621.1,  Hypersonic  Zero  Pressure  Gradient 
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readings,  thereby  introducing  the  least  possible  instrument  error.  Skin  friction 
values  were  measured  using  a  floating  element  shear  stress  gage.  As  can  be 
seen  in  figure  16,  Mach  number  is  underpredicted  by  about  15  percent  by  all  the 
theories  when  this  measured  shear  stress  level  is  input  to  the  WALAW  program. 

A  value  of  tw  approximate.lv  20  percent  higher  would  give  much  better  agreement, 
as  shown  in  the  figure.  Similar  results  are  obtained  for  the  static  tempera¬ 
ture  prediction,  as  shown  in  Figure  17.  Static  temperature  is  overpredicted 
using  the  reported  wall  shear;  a  tw  value  20  percent  larger  gives  more  corn  't 
predictions  with  all  theories  considered.  More  will  be  said  of  the  agreement 
with  Lee's  data  later  in  this  report. 

Lee's  higher  stagnation  pressure  experiments  show  better  agreement  with 
theoretical  predictions  using  the  reported  shear  values.  Figure  18  shows  Mach 
number  data  and  predictions  for  a  P  ■  150  psia  run,  the  highest  pressure 
reported  by  Lee.  Agreement  between  Mach  number  data  and  profiles  is  excel¬ 
lent  for  all  three  theories,  as  are  the  static  temperature  predictions 
(Figure  19) .  The  effect  of  a  20  percent  change  in  wall  heat  flux  is  also 
shown  in  this  figure,  in  order  to  establish  the  sensitivity  of  the  one-dimen¬ 
sional  solution  to  this  input  quantity.  Thus,  it  can  be  concluded  that  for 
incompressible  flows  or  compressible  flows  up  to  a  Mach  number  of  5,  with 
or  without  heat  tr  n-.fer,  all  three  turbulent  models  do  a  satisfactory,  and 
a  very  similar,  job  of  predicting  profiles  in  the  near  wall  region. 

As  discussed  earlier,  Aerotherm's  mixing  length  law  has  been  verified 
by  comparisons  with  Kendall's  (refererce  66)  low-speed  blown  boundary  layer 
data  in  another  report  (reference  2) .  These  comparisons  were  rerun  with 
the  WALAW  program  to  confirm  their  validity  and  also  to  obtain  simultaneous 
predictions  for  the  Cebeci  and  the  Bushnell-Beckwith  models.  A  small  error 
was  found  in  the  original  (reference  2)  comparisons;  however,  the  results 
remain  essentially  the  same,  as  can  be  seen  in  Figures  20,  21,  and  22.  All 
results  were  obtained  using  the  wall  shear  values  reported  by  Kendall,  deter¬ 
mined  by  a  wall  pressure  profile  technique  (as  opposed  to  a  d6/dx  technique). 

At  the  mooaat  blowing  rate  of  F  ■  0.002  (Figure  20),  agreement  between  the 
Aerotherm  and  Cebeci  models  and  the  experimental  data  is  excellent.  The 
Bushnell-Beckwith  model  appears  to  react  too  strongly  to  the  blowing  rate,  as 
evidenced  by  the  very  small  extent  of  the  laminar  sublayer.  A  higher  input 
value  would  shift  all  predictions  upward,  however  (see  Figure  16).  In 
i-igure  21,  results  for  a  stronger  blowing  rate  of  F  «  0.003  are  shown.  At  this 
blowing  rate,  the  Aerotherm  model  shows  the  leant  effect  of  blowing,  while 
the  Pushnell  model  shows  the  most.  Agreement  is  good  for  Model  A  and 
Model  C.  Results  for  strong  blowing  are  shown  in  Figure  21.  Using 
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Linear-Log  Velocity  Ratio  Profile 

MIT,  Kendall,  Run  C-2xl0  s  -50,  Low-Speed  Flow  With  Blowing 


Figure  21.  Linear- Log  Velocity  Ratio  Profile 

MIT,  Kendall  Run  C  -  3  x  10“3  -  50,  Low-Speed  Flow  With  Blowing 
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the  reported  wall  shear  stress  and  blowing  rate,  orly  the  Aerotherm  model 
predicts  the  velocity  profile  accurately.  The  Cebeci  and  Bushnell  models 
result  in  nearly  idencical  profiles,  giving  an  error  of  50  percent  in  the 
law  of  the  wall  region.  The  uncertainty  in  wall  shear  can  more  than  account 
for  this  entire  error,  however. 

The  Stanford  experimental  data  (reference  33)  provides  an  independent 
check  on  the  low  speed  blowing  results  presented  above,  since  the  test  con¬ 
ditions  are  nearly  identical.  Froestream  velocities  of  approximately  50  ft/sec 
with  no  pressure  gradient  were  used  for  the  reference  33  experiments,  with 
blowing  rates  of  F  =  0.0  to  0.008.  Figure  23  illustrates  th3  good  agreement 
obtained  between  experiment  and  theory  for  a  no-blowing  case.  Figures  24,  25, 
and  26  present  data  and  predictions  for  nominal  blowing  rates  of  F  =  0.002, 
0.004,  and  0.008,  respectively.  Tne  theory  trends  are  exactly  as  shown  in  the 
Kendall  data  comparisons.  Models  A  and  C  show  nearly  the  same  results  at  low 
blowing,  but  diverge  at  stronger  blowing  rates.  Models  B  and  C  give  the  same 
results  with  strong  blowing.  Using  the  reported  wall  shear  values,  the  data  in 
general  follow  the  trends  predicted  by  the  Bushnell-Beckwith  model  more  closely 
than  the  Cebeci  model.  Both  of  these  models  do  a  better  job  than  the  Aerotherm 
model.  Thus,  any  conclusions  which  may  have  been  drawn  from  the  Kendall  data 
comparisons  must  be  questioned  in  the  light  of  the  Simpson  comparisons.  The 

difficulty  lies  in  the  requirement  of  an  accurate  C-  value,  which  is  uncertain 

*  1 

to  +  100  percent  or  more  for  strongly  blown  flows. 

One-dimensional  theory  and  data  comparisons  were  also  carried  out  for 
the  supersonic  blown  boundary  layer  experiments  of  Jeromin  (reference  48) . 
Jeromin  experienced  difficulties  with  axial  pressure  gradients  and  non-two- 
dimensionality  of  the  flow,  therefore  his  data  must  be  considered  less  depend¬ 
able  than  the  subsonic  cases.  Drag  coefficients  were  inferred  from  the  integral 
momentum  equation,  from  Stevenson's  wall  law  technique  (reference  81),  and  by 
a  transformation  method.  In  the  predictions  for  several  of  Jeromin' s  runs  shown 
in  Figures  27  to  30,  a  single  tw  value  was  used  if  the  three  Cf  determination 
methods  were  in  agreement.  More  than  one  xw  value  was  used  (and  is  illustra¬ 
ted  in  the  figure)  if  the  methods  gave  widely  divergent  results.  Figure  27 
shows  the  type  of  agreement  obtained  for  an  unblown,  Mach  2.5  flow.  Agreement 
with  the  predictions  is  fairly  good  at  the  lower  indicated  shear  stress  value, 
although  not  as  good  as  the  Coles'  comparisons  in  supersonic  flow  (Figure  15). 
Moderate-to-strong  blowing  at  this  Mach  number  (Figure  28)  gives  poor  agreement 


*  The  Kendall  wall  pressure  profile  technique  for  determining  Cf  was  applied  to 
several  of  Simpson's  runs,  but  no  identifiable  trends  in  the  alternate  values 
for  Cf  were  found. 
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Figure  23.  Linear-Log  Velocity  Ratio  Profile 

Stanford,  Simpson  Run  2367-3,  Low-Speed  Flow  Without  Blowing 


Figure  24.  Linear-Log  Velocity  Ratio  Profile 

Stanford,  Simpson  Run  122366-2,  Low-Speed  Flow  With  Blowing 


Ratio  Profile 

in  12206C-3,  Low-Spoed  Flow  With  Blowing 


Figure  26.  Linear-Log  Velocity  Katio  Profile 

Stanford,  Simpson  Run  121966-4,  Low-Speed  Flow  With  Blowing 


.40,  Supersonic  Flow  without  ilowing 


Supersonic  Flow  With  Blowing 
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with  Model  A  and  good  agreement  with  Models  B  and  C.  Results  for  a  Mach  3.5 
flow  are  shown  in  Figures  29  and  30.  With  no  blowing,  agreement  is  again  fairly 
good  for  all  theories.  A  very  strong  blowing  case  was  selected  for  the  last 
comparison,  Figure  30.  As  with  the  Simpson  strong  blowing  case,  agreement  with 
Models  B  and  C  is  good,  while  Model  A  gives  poor  agreement.  However,  an  un¬ 
certainty  in  tw  of  100  percent  for  this  case  is  not  unlikely. 

The  preceding  comparisons  have  shown  that  three  leading  turbulent  models 
for  near  wall  calculations  are  in  substantial  agreement  for  unblown  and  small- 
to-modeiately  blown  flows.  At  larger  blowing  rates,  with  the  wall  shear  level 
specified,  the  Aerotherm  model  is  in  substantial  disagreement  with  the  other 
two  models  and  with  the  experimental  data  of  Simpson  a  cf  Jeromin.  On  the 
other  hand,  the  Aerotherm  model  is  in  agreement  with  the  strong  blowing  ex¬ 
perimental  data  of  Kendall,  which  is  not  matched  by  the  Cebeci  and  Bushr.ell 
•theories.  Thus,  the  choice  of  a  "best”  turbulent  model  is  not  at  all  clear. 

The  question  of  which  model  is  best  cannot  actually  be  answered  with 
these  one-dimensional  theory  comparisons,  since  the  predictions  are  so  completely 
dependent  on  the  wall  shear  rate.  For  blown  boundary  layers,  the  wall  shear  is 
generally  highly  uncertain.  Thus,  insufficient  information  is  available  to 
make  a  choice  between  the  three  models  on  purely  technical  comparison  grounds. 

The  basis  for  the  choice  which  was  made  is  presented  in  the  next  subsection. 

3.  THE  SELECTED  TURBULENT  MODEL 

a.  Wall  Law 

The  previous  comparisons  have  centered  on  the  wall  law  region,  since 
this  area  of  turbulence  modeling  is  perhaps  the  most  importan'  and  includes  the 
widest  disagreement  between  theories.  The  comparisons  of  Section  III. 2  have 
shown  substantial  agreement  between  theories  for  flows  with  no  blowing  or  "small" 
blowing,  but  wide  disagreement  for  flows  with  "strong"  blowing. 

In  order  to  place  the  "small"  blowing  and  "strong"  blowing  terminology 
in  perspective,  a  practical  example  is  presented  here.  Consider  a  10°  half- 
angle  cone  flying  at  sea  level  at  Mach  10.  Assume  an  ablating  carbon  or  phe¬ 
nolic  carbon  surface  on  the  cone,  such  that  the  wall  temperature  is  about  6000°R. 
With  all  turbulent  flow,  the  "no  blowing"  drag  coefficients  at  the  1-fcot  and 
5-foot  body  stations  are  approximately 

Cf  ■  0.00042  §  x  »  1  foot 

Cj  ■  0.00C33  @  x  «■  5  feet. 
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These  values  were  calculated  using  the  Schultz-Grunow  formula  corrected  for 
compressibility j 


0.37 


i  (l°g10Rex) 


2.51 


(31) 


(32) 


*  1  +  0 
e 


(33) 


For  a  carbon  nose  tip  or  heat  shield,  the  ablation  rate  is  characterized  by 
plateau-like  behavior  at  a  B’  =  m  /(p  u C„)  =0.2  (reference  82).  Thus,  assuming 

"  C  6  6  n 


C„  -  Cf/2, 


/C£ 

7p5TT  =  °*2  4 


(34) 


which  results  in  a  very  small  blowing  rate.  A  similar  behavior  is  obtained 
with  phenolic  carbon  at  a  of  approximately  0.6.  Thus,  the  blowing  rates  tor 
typical  heat  shield  materials  are  quite  small  under  ordinary  flight  conditions. 

In  light  of  the  fact  that  there  are  only  small  differences  between  modals 
at  zero  or  small  blowing,  and  that  small  blowing  conditions  are  typical  for 
flight  cases,  the  logical  choice  of  a  wall  region  turbulent  model  for  BLIMP  com¬ 
parisons  was  the  reference  1  model.  The  reference  1  model  is  already  incor¬ 
porated  in  the  BLIMP  code, which  minimized  the  amount  of  programming  time  re¬ 
quired.  The  constants  within  this  model  were  held  at  the  incompressible  flow 
values  of 


k  -  0.44 
m 


yZ  •  11.823 

U 


(35) 

(36) 


The  effects  of  changes  in  these  constants  will  be  demonstrate*  in  Section  IV. 
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b.  Wake  Law 

There  has  been  very  little  discussion  in  the  literature  about  the 
advantages  or  disadvantages  of  a  Clruser  type  wake  law,  as  used  in  Models  A  and 
C  compared  to  a  boundary  layer  thickness  dependent  mixing  length  law,  as  used 
in  Model  B.  All  that  seems  to  be  required  in  the  outer  portions  of  equilibrium 
type  flows  is  that  a  roughly  correct  eddy  viscosity  number  be  arrived  at.  This 
is  borne  out  by  the  BLIMT  results  shown  in  figure  3.1  ,  where  the  constant  in  the 
expression 


em  =  (const. ) ug5£  (37) 

was  set  at  both  0.016  and  0.018,  all  other  features  of  the  reference  1  model 
remaining  the  same.  This  12.5  percent  change  in  the  outer  boundary  layer  tur¬ 
bulence  coefficient  over  a  distance  of  2.75  feet  in  running  length  produced 
virtually  no  change  in  the  velocity  profile.  The  Clauser  type  formulation  in 
terms  of  the  kinematic  or  velocity  displacement  thickness  6|,  seems  to  be 
adequate  for  equilibrium  compressible  flows,  therefore  the  reference  1  wake 
model  was  adopted  for  the  final  data  coifiparisons  using  the  0.018  constant. 

c.  Turbulent  Prandtl  and  Schmidt  Number 

The  turbulent  Prandtl  and  Schmidt  numbers  relate  the  turbulent  trans¬ 
port  of  energy  and  chemical  species  to  the  turbulent  transport  of  momentum. 

While  there  is  only  a  weak  basis  for  doing  so,  it  is  common  practice  to  assume 
a  constant  jatio  between  the  transport  coefficients.  Experimental  data  on 
actual  turbulent  Prandtl  numbers  and  their  variations  through  the  boundary  layer 
are  becoming  available  (references  33,  84,  and  85).  The  evidence  indicates  that 
Prfc  should  vary  from  ~  0.5  in  the  outer  portions  of  the  boundary  layer,  to  near 
unity  at  y/<$  =  0.1,  to  ~  1*5  or  2.0  near  the  laminar  sublayer.  Results  are 
primarily  for  air  boundary  layers.  While  these  results  are  certainly  relevant 
to  the  BLIMP  turbulent  model,  the  development  and  incorporation  of  a  Prfc  model 
which  has  reasonable  validity  for  multicomponent  chemically  reacting  flows  was 
beyond  the  scope  of  the  current  effort.  Thus,  a  constant  Prt  model  was  used 
for  the  current  studies.  A  similar  situation  exists  for  turbulent  Schmidt 
number,  although  fewer  data  are  available  (e.g.,  reference  86). 

Numerical  studies  with  the  BLIMP  code  and  the  reference  1  turbulent  model 
have  shown  that  a  turbulent  Prandtl  number  of  0.9  results  in  the  generally 
accepted  value  of  0.88  for  the  flat  plate  recovery  factor  in  air.  These  results, 
shown  in  Figure  32,  were  generated  with  BLIMP  by  demanding  zero  heat  flux  to 
the  wall  as  a  solution  boundary  condition.  Thus,  a  constant  value  of  Prfc  ■  O.y 
was  used  for  boundary  layer  predictions  with  the  current  turbulent  model. 


65 


0.10 


0.09 


0.08 


0.07 


0.06 


I  I  I 


TYPICAL  HYPERSONIC  PROFILE 


predict: 

M  ■*  6.0  _l 

EON 

r _ i 

flat  plate 

X  -  2.75  ft 

*'  I 

1 

■ 

'blimp  PREDICTION  WITH 
*  0.016 


-  BLIMP  PREDICTION  WITH 

-  e  ■  0.018  u,it 

m  1  1 


0.05 


Figure  32.  Effect  of  Turbulent  Prandtl  Number  on  Recovery  Factor  at 


AFWL-TR-71-57 


SECTION  IV 

PREDICTIONS  AND  COMPARISONS  WITH  DATA 

1.  PRESENTATION  OF  COMPARISONS 

The  comparisons  between  BLIMP  predictions  and  measured  data  for  the  17 
separate  data  sets  listed  in  Table  V  are  presented  in  this  section.  The 
results  are  listed  in  the  order  of  Table  V.  They  consist  of  (1)  a  brief  des¬ 
cription  of  the  BLIMP  input*  as  related  to  the  data  and  (2)  graphical  results, 
where  applicable,  in  the  form  of  linear  Mach  number  ratio,  total  temperature 
ratio,  and  velocity  ratio  profiles  as  functions  of  y,  velocity  ratio  profiler 
as  functions  of  y  in  linear-log  coordinates,  and  linear  plots  of  momentum  thick¬ 
ness  Reynolds  number,  Re0,  and  skin  friction  coefficient,  Cf/2,  versus  stream- 
wise  dimension,  x.  Symbolism  has  been  standardized  for  all  plots.  Measured 
data  are  represented  by  circled  points;  where  more  than  one  set  exists  on  a 
single  plot,  symbols  are  noted  on  that  plot,  BLIMP  predictions  started  far 
upstream  of  the  region  of  interest  using  an  approximate  starting  profile  (a 
"zero  start")  are  presented  as  solid  curves;  predictions  started  at  the  first 
data  station  using  an  actual  data  profile  (a  "data  start")  are  presented  as 
dashed  curves.  Any  deviations,  modifications,  or  additions  to  this  format  are 
noted  on  the  figure.  In  general,  only  one  linear-log  velocity  plot  is  shown, 
usually  the  next  to  last  of  those  plotted  in  linear  coordinates. 

a.  Coles 1  Comparisons 

(1)  Comments  on  BLIMP  Input 

Coles'  data  consist  of  a  single  profile  with  measurements  at 
three  other  axial  stations.  Freestream  and  wall  conditions  were  known  to  be 
effectively  constant,  thus  BLIMP  was  started  from  the  leading  edge  of  the  flat 
plate.  The  use  of  a  fence  trip,  however,  introduces  an  unknown  initial  thick¬ 
ness  which  is  equivalent  to  the  virtual  origin  being  upstream  of  the  leading 
edge.  Since  the  data  start  option  is  not  applicable  with  only  a  single  station 
of  profile  data  available,  the  only  means  of  accounting  for  the  trip  effect  was 
to  match  the  profiles.  The  momentum  thickness  was  selected  as  the  appropriate 


*  A  more  complete  discussion  of  computer  code  setup  for  these  problems,  includ¬ 
ing  zero  starts  and  data  startB,  is  included  in  the  appendix. 
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matching  parameter.  A  BLIMP  profile  at  the  desired  value  of  0  was  obtained 
by  first  running  the  program  with  sufficient  stations  to  bracket  the  data  value, 
then  rerunning  the  program  including  the  interpolated  value  of  x  as  one  of  the 
output  streamwise  stations  (see  Figure  39).  Thus  the  profile  comparisons  of 
Figures  33-38  compare  data  taken  at  x  *  1.69  feet  to  predictions  at  2.46  and 
2.18  feet  for  M  ■  2.5  and  3.6,  respectively. 

(2)  Comments  on  Results 

Mach  number  and  velocity  ratio  comparisons  in  Figures  33-38  for 

both  runs  are  relatively  good.  Temperature  profiles  are  not  included  since 

temperatures  were  not  measured.  In  general,  the  prediction  indicates  higher 

values  of  Mach  number  and  velocity  in  the  mid-ranges  of  the  profile  and  lower 

values  beyond  u/ug  »  0.8  to  0.9.  There  is  a  slight  indication  of  an  inflection 

point  in  the  Mach  number  data  profile  which  the  prediction  does  not  seem  to 

include.  Maximum  absolute  differences  in  ratio  are  0.04  to  0.05.  The  plot 

of  Re0  in  Figure  33  is  shown  only  to  indicate  the  matching  procedure.  Finally 

the  C^/2  plot  in  Figure  40  shows  the  BLIMP  predictions  and  data  both  referenced 

from  the  beginning  of  the  plate.  C^/2  predictions  for  the  6-matched  stations 

are  7.5  and  8.5  percent  lower  than  the  data  for  the  Mach  2.5  and  3.6  cases, 

respectively.  It  should  be  noted,  however,  that  improvements  in  the  nodal 

distribution  to  be  discussed  in  Section  IV. 2. b  account  for  most  or  all  of 

this  error.  Nodal  spacing  was  discovered  to  be  particularly  important  for 

dQ 

8 -matched  cases,  in  that  a  small  improvement  in  brings  about  a  new  match 
point,  which  in  turn  gives  a  significantly  different  C^. 

b.  NOL  Comparisons 

NOL  data  sets  are  comprised  of  profiles  measured  at  four  stations 
(Lee)  and  five  stations  (Brott)  with  skin  friction  measurements  several  inches 
upstream  of  each  profile  station.  Measurements  were  taken  over  a  streamwise 
interval  from  approximately  4  to  7,5  feet  measured  from  the  nozzle  throat: . 

Each  profile  was  measured  during  a  separate  test  run,  thus  stagnation  conditions 
varied  slightly  (less  than  one  percent)  from  profile  to  profile.  Comparison 
to  a  single  BLIMP  run,  which  is  based  on  constant  stagnation  conditions,  is 
subject  to  some  small  error  for  this  reason.  In  general,  the  BLIMP  input  of 
freestream  conditions  was  based  on  the  measured  Mach  number  and  average 
stagnation  conditions  for  the  four  or  five  reported  profiles.  Wall  temperatures 
were  taken  as  the  reported  values.  As  a  consequence  of  this  approach,  the 
pressure  gradients  at  a  given  test  station  are  slightly  different  for  the  predic¬ 
tion  than  actually  existed.  Axial  pressure  gradients  are  not  discussed  by 
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Figure  37.  Velocity  Ratio  Profile 

Coles  Run  #20,  Supersonic,  Zero  Pressure  Gradient  FI 
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Figure  36.  Linear-Log  Velocity  Ratio  Profile 

Coles  Run  #20,  Supersonic,  Zero  Pressure  Gradient  Flow 
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Figure  39.  Momentum  Th.-okncss  Reynolds  Number  vs  Streamwise  Location 
Celes  Runs  at  luich  2.6  find  3.7 
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Figure  40.  Skin  Friction  Coeff  cient  vs  Streamwise  Location 
Coles  Runs  at  Mach  2.6  and  3.7 
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Lee,  et  al . ,  except  to  note  that  the  48-inch  station  was  located  in  the 
acceleration  region,  55  inches  being  the  beginning  of  the  constant  freestr  Rm 
flow  where  Mach  number  variations  are  constant  within  ±  0.75  percent.  There 
is  no  indication  of  local  gradient  conditions  at  either  the  skin  friction 
gages  or  profile  station.  Srott  does  list  such  data;  however,  there  is 
currently  no  way  to  model  the  pressure  gradients  directly  in  the  BLIMP  code. 
Pressure  gradients  at  the  prediction  stations  are  determined  by  a  quadratic 
curve  fit  of  the  input  pressure  distribution  which  depends,  of  course,  on  the 
interval*  between  stations  as  well  as  the  pressures. 

(1)  Data  of  Lee,  et  al. 

(a)  Comments  on  BLIMP  Input 

To  show  the  effects  of  pressure  gradients,  the  5  atmosphere 
run  of  Lee  was  set  up  as  described  above,  while  constant  freestream  conditions 
equal  to  the  average  of  those  reported  at  the  four  stations  were  used  for  the 

10  atmosphere  run.  Both  cases  were  run  with  the  data  start  option  since  con- 

* 

ditions  upstream  of  the  first  profile  station  were  not  reported.  Several 
"approximate"  zero  start  runs  were  attempted  for  the  10  atmosphere  run,  one 
using  the  nozzle  profile  reported  by  Brott  and  the  other  assuming  constant  con¬ 
ditions  throughout.  These  results  bracketed  the  measured  values  of  6  only  in 
a  gross  sense  (see  Figure  51).  However,  as  shown  in  Figure  52,  these  approxi¬ 
mate  zero  start  runs  do  define  limits  on  the  value  of  Cf/2  within  which  a  pre¬ 
diction  based  on  the  actual  upstream  conditions  can  be  expected  to  fall. 

(b)  Comments  on  Results 

Figures  41-46  and  47-51  show  the  results  for  the  5  and  1C 
atmosphere  runs,  respectively.  In  both  cases,  the  predicted  and  measured  profiles 
for  Mach  number,  velocity  and  temperature  diverge  from  the  data  as  the  flow  pro¬ 
ceeds  down  the  plate  (after  the  matched  data  start  beginning) .  The  predictions 
exceed  the  data  in  the  mid-region  of  the  Mach  number  and  velocity  profiles  by 
up  to  30  percent.  In  Figures  41  and  47  the  Mach  number  data  profiles  appear  to 
undergo  a  significant  change  in  shape  that  is  not  reflected  by  the  predictions. 

In  the  Re0  plots  of  Figures  45  and  51,  the  predicted  growth,  starting 
from  the  matched  value  at  the  first  station,  is  less  than  that  measured.  This 
is  consistent  with  the  profile  plots  which  show  progressively  "thicker" 


*  Private  communication  with  NOL  personnel  indicated  these  data  were  not 
measured. 
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Figure  44.  Linear-Log  Velocity  Ratio  Profiles 

SiOL ,  Lee,  et  al.,  5  Atm,  Hypersonic,  Zero  Pressure  Gradient  Flow 


Figure  *5,  Momentum  Thickness  Reynolds  Number  vs  Streamwise  Location 

KOL,  Lee,  et  al.,  5  Atm,  Hypersonic,  Zero  Pressure  Gradient  Flow 


FIGURE  47.  !WCH  IM1BER  PATlO  PROFILES 

HOt.  LEE.  ET  AL..  10  ATR.  HYPERSOIJIC,  ZERO  PRESSURE  GRADIENT  FLOW 


TBIPtRATURE  RATIO  PROFILES 
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Figure  50.  Linear-Log  Velocity  Ratio  Profiles 

NOL,  Lee,  et  al.,  10  Atm.,  Hypersonic,  Zero  Pressure  Gradient  Flow 


(ASSUMED  dP /'X  VARIA- 


Figure  51.  Momentum  Thickness  Reynolds  Number  vs  Streamwise  Location 

NOL,  Lee,  et  al.f 10  Atm,  Hypersonic,  Zero  Pressure  Gradient  Flow 


(ASSUMED  dP/dx  VARIA 


Figure  52.  Skin  Friction  Coefficient  vs  Streamwise  Location 

NOL,  Lee,  et  al.,  10  Atm,  Hypersonic,  Zero  Pressure  Gradient  Flow 
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measured  boundary  layers  in  terms  of  0  even  though  the  predicted  and  measured 
6  thicknesses  remain  essentially  equal  (see  Figures  43  and  49).  The  consis¬ 
tency  in  the  data  appears  doubtful  when  one  considers  the  C^/2  plots  of 
Figures  46  and  52.  For  both  runs,  the  C^/2  predictions  exceed  the  measured 
values  after  the  first  station  (which  is  still  within  the  accelerating  region 
of  the  nozzle)  by  16  to  24  percent.  This  is  directly  opposite  to  the  trend 
seen  in  the  Re@  plots  (Figures  45  and  51).  Considering  the  momentum  integral 
equation,  which  for  zero  pressure  gradient  and  no  blowing  reduces  to 


%  «  p.  (38) 

2  3x 

J  A 

it  is  clear  that  if  C^/2  is  overpredicted,  should  also  be  overpredicted. 

The  data,  of  course,  result  from  skin  friction  gages,  rather  than  a  solution 
of  the  momentum  integral  equation.  Calculation  of  Cf/2  using  equation  (38) 
and  measured  values  of  0  and  x  results  in  values  between  8  and  9  x  ltf  as  com¬ 
pared  to  skin  friction  data  between  3.6  and  4.0  x  10*'  in  the  uniform  flow 
region.  On  the  other  hand,  the  BLIMP  predicted  Cf/2  and  d0/dx  are  equal. 

In  view  of  the  above  comparisons,  the  Lee  dal  a  appear  somewhat  question¬ 
able.  This  is  not  to  say  the  predictions  are  100%  correct,  but  only  that 
internal  inconsistencies  in  the  data  such  as  three-dimensional  flow  are  sus¬ 
pected.  Further  comment  will  be  made  in  Brott  comparisons  to  follow. 

The  effect  of  the  "input"  pressure  gradients  in  the  Lee  data  is  most 
apparent  in  the  Re0  and  Cf/2  plots  of  Figure  45,  46,  51,  and  52.  In  the  first 
set  (5  atm,  Figures  45,  46)  pressure  variations  are  considered  in  the  predict¬ 
ion  while  in  the  latter  (Figures  51,  52)  they  are  not.  The  axial  pressure 
gradient  distribution  computed  by  BLIMP  is  negative  between  3.77  and  4.7  feet, 
positive  up  to  7.0  feet  and  egative  thereafter.  The  effect  on  Re0  is  seen 

in  Figure  45;  d9/dx  la’  in  regions  of  negative  pressure  gradient  and 

* 

smaller  in  the  positive  gr«_  *ent  regions.  The  momentum  integral  equation 
accounting  for  pressure  gradients  has  the  form 


C 


f 


~T 


(39) 


*  This  behavior  in  Cf/2  is  directly  opposite  that  for  subsonic  flow  due  to 
the  change  of  sign  of  the  bracket  term  in  equation  39  at  high  Mach  numbers. 
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Table  /III  lists  the  values  of  each  of  these  terms  at  each  of  seven  BLIMP 
output  stations. 

TABLE  VIII 

MAGNITUDES  OF  MOMENTUM  INTEGRAL  EQUATION  TERMS 


Pressure 

Term 

Cf 

T 

3.77 

13.00  *  lo"4 

8.25  *  10"4 

4.75  *  10'4 

3.82 

19.18 

14.00 

5.18 

3.90 

18.58 

13.60 

4.98 

4.20 

10.24 

5.45 

4.79 

4.77 

4.41 

-0.32 

4.73 

5.77 

0.04 

-4.61 

4.65 

7.60 

9.56 

4.80 

4.76 

From  Table  VIII  and  Figure  46 ,  it  is  apparent  that  the  predicted  value 
of  Cf/2  experiences  only  minor  fluctuations  even  though  the  value  of  the  pres¬ 
sure  term  is  oscillating  with  an  amplitude  several  times  the  value  of  Cf/2. 

That  is,  the  input  pressure  gradient  is  a  very  important  term  in  the  momentum 
integral  equation;  however,  adjustments  to  account  for  its  large  variations 
occur  in  the  e  growth  rate  (d6/dx)  rather  than  in  Cf/2.  The  initial  spike  in 
Cf/2  is  due  to  the  incorrect  input  pressure  gradient.  The  spike  has  no  real 
significance  since  the  second  through  the  fourth  stations  were  added  only  for 
purposes  of  the  data  start  option,  and  the  hand-interpolated  pressure  values 
resulted  ir.  locally  high  gradients.  These  pressure  gradients  offer  no  obstacle 
to  the  computer  solution,  and  have  no  lasting  effect  on  the  downstream  solutions; 
therefore,  the  solution  was  not  rerun.  The  results  do  emphasize  the  experimental 
difficulties  associated  with  accurate  evaluation  of  local  values  through  use 
of  the  momentum  integral  (equation 

In  Figure  51  the  momentum  growth  is  essentially  linear  as  expected  for 
a  constant  axial  pressure  and  nearly  constant  throughout  the  region.  The 
variation  in  Cf/2  in  Figure  52  is  due  to  several  factors.  The  initial  rise 
is  the  rapid  recovery  of  the  profile  next  to  the  wall  following  the  data  start 
profile  as  described  in  the  appendix.  The  slow  variation  up  to  x  »  6.0  feet  is 
believed  to  be  due  to  the  readjustment  of  the  outer  portion  cr  the  profile 
which  does  have  a  small  effect  on  wall  gradients.  Downstream  of  x  *  6.f  ;'eet, 
this  readjustment  appears  to  be  complete  and  C^/2  decreases  slowly  an  expected 
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for  a  uniform  flow.  These  observations  ere  further  substantiated  by  comparison 
to  the  two  solid  curves  which  indicate  the  approximate  zero  start  BLIMP  predic¬ 
tions,  assuming  the  upstream  flow  conditions  noted. 

(2)  Data  of  Brott,  et  al. 

(a)  Comments  on  BLIMP  Input 

Vhe  Mach  number  distribution  from  the  throat  to  the  first 
profile  station  is  presented  in  reference  6  for  the  design  and  the  measured  val¬ 
ues.  Using  this  information,  zero  start  BLIMP  predictions  were  made  assigning 
x  =  0  at  the  nozzle  throat.  For  additional  comparison,  "data  start"  runs  start¬ 
ing  from  the  first  reported  profiles  at  x  =  47  inches  were  made.  As  witn  the  5 
atmosphere  Lee  run  above,  freestream  pressures  were  assigned  to  equate  the 
local  measured  and  input  freestream  to  stagnation  pressure  ratios. 

(b)  Comments  on  Results 

Figures  53-58  and  59-64  contain  the  results  of  the  5  and  10 
atmosphere  stagnation  pressure  runs,  respectively.  Results  are  plotted  for 
the  zero  start  and  data  start  predictions  described  above.  Profiles  fron  the 
data  start  prediction  for  both  runs  tend  to  readjust  in  shape  moving  down  the 
plate.  By  the  final  station,  which  is  nearly  3  feet  downstream  from  the  first, 
these  profiles  are  quite  similar  to  the  zero  start  profiles.  This  then  is 
a  measure  of  the  duration  of  the  recovery  process.  These  figures  show  that 
the  adjustment  of  the  Mach  number  occurs  more  rapidly  than  that  of  total  temp¬ 
erature  and  velocity. 

In  comparison  with  the  data,  both  forms  of  the  prediction  tend  to  (1) 
indicate  higher  values  of  Mach  number,  total  temperature,  and  velocity  in  the 
mid-range  of  the  profiles  up  to  Mach  number  ratios  of  0.8  and  velocity  ratios 
of  0.9  and  (2)  underpredict  in  the  far  wake  region,  particularly  the  Mac.n 
number.  As  with  Lee  above,  the  boundary  layer  thicknesses,  6,  are  essentially 
equal. 

The  momentum  thickness  comparisons  of  Figures  57  and  63  are  quite  good 
for  both  starting  conditions.  Small  differences  in  0  at  the  first  data  point 
are  apparently  due  to  differences  between  the  Aerotherm  and  NOL  curvefit  of 
the  profile  data.  This  good  overall  agreement  is  in  marked  contrast  to  the 
results  of  the  Lee  comparisons.  Figures  58  and  64  show  excellent  agreement 
between  predictions  and  data  for  Cf/2.  Note  that  there  is  no  sudden  change  in 
Cf/2  after  the  first  station  for  the  data  start  predictions,  indicating  that 
the  input,  i.e.,  the  measured  wall  gadient,  was  consistent  with  the  BLJ'IP  solu¬ 
tion  for  the  local  conditions.  The  more  gradual  dips  shown  are  probably  due  to 
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Figure  56.  Linear-Log  Velocity  Ratio  Profiles 

NOL,  Brott,  et  al.,  5  Atm,  Hypersonic,  Negative  Pressure  Gradient  Flow 
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Figure  62.  Linear-Log  Velocity  Ratio  Profiles 

NOr,,  Brott,  et  al.,  10  Atm,  Hypersonic,  Negative  Pressure  Gradient  Flow 
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Figure  64.  Skin  Friction  Coefficient  vs  Streamwise  Location 

NOL,  Brott,  et  al.,  10  Atm,  Hypersonic,  Negative  Pressure  Gradient  Flow 
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the  overall  profile  adjustnvent  taking  place  but  may  also  be  a  result  of  the 
modeling  of  the  pressure  gradient  in  this  region. 

The  log  velocity  plots  of  Figures  56  and  62  also  show  good  agreement 
over  the  law  of  the  wall  region.  The  two  data  points  nearest  the  wall  in  both  of 
these  'igures  are  assumed  to  be  in  error  since  they  are  inconsistent  with  the 
skin  friction  balance  data  and  the  BLIMP  preaiction.  These  data  points  were 
taken  with  a  very  small  probe  (0.005  inch  high)  at  very  low  pressures,  ^nd 
are  subject  to  numerous  sources  of  error.  Referring  back  to  the  Lee  log 
velocity  plots  (Figures  44  and  50),  near  wall  measurements  show  better  agree¬ 
ment  with  theory,  perhaps  since  a  larger  (0.016  inch)  probe  was  used.  When 
the  near  wall  readings  from  three  stations  are  superimposed  (Figures  44  and 
50),  the  5  atmosphere  case  shows  good  agreement  while  the  10  atmosphere  case 
does  not.  This  merely  confirms  the  idea  that  the  near  wall  data,  even  with 
the  larger  probe,  are  questionable.  Farther  away  from  the  wall,  over  the  bulk 
of  the  profile,  the  Lee  data  are  in  much  poorer  agreement  with  theory  than 
Brott's  data. 

In  conclusion,  the  Brott  data  appear  to  be  much  more  consistent  in  terms 
of  profile  shape,  0  growth,  and  reported  C^/2  values,  both  internally  and  in 
comparison  to  the  predictions.  The  Lee  data  lack  this  internal  consistency, 
particularly  when  d6/dx  is  compared  to  the  Cf/2  values  Also  in  the  Lee 
experiments,  the  rapidly  changing  profile  shape  in  an  essentially  ’’similar" 
flow  region  is  disturbing. 

c.  TRW  Comparisons 

(1)  Comments  on  BLIMP  Input 

Since  tunnel  operation  was  continuous,  variations  in  stagnation 
and  freestream  conditions  from  profila  to  profile  were  quite  small.  Constant 
input  conditions  were  determined  as  averages  of  the  reported  station  values. 

Wall  temperatures  were  input  as  measured. 

An  experimental  difficulty  that  might  have  affected  the  profile  data  was 
the  formation  of  a  frost  layer  on  the  cold  wall  section.  This  was  not  detected 
until  after  the  step-up  run.  h  change  was  made  in  procedure,  namely,  data  were 
taken  during  only  the  first  15  to  20  minutes  once  the  wall  was  cooled.  The 
model  was  then  warmed  and  recooled  until  measurements  were  complete.  This  pro¬ 
cedure  was  employed  for  the  step-down  tests.  Due  to  the  frost  formation,  Stan¬ 
ton  tube  data  in  cold  wall  regions  were  judged  to  be  erroneous  and  thus,  not 
reported. 


*  Personal  communication  with  Roland  Lee  at  NOL  indicated  that  correct! >n  curves 
for  near  wall  measurements  like  these  are  currently  being  generated. 
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Figures  65  and  71  present  the  axial  wall  temperature  distributions  for 
the  step-up  and  step-down  runs.  BLIMP  comparisons  were  made  at  the  stations 
marked  by  checks.  For  the  zero  start  predictions,  freestream  and  wall  condi¬ 
tions  upstream  of  the  first  reported  station  were  assumed  constant  and  equal  to 
the  values  at  the  first  station.  In  handling  step  changes  in  streamwise  prop¬ 
erties,  BLIMP  contains  an  option  to  treat  streamwise  derivatives  as  two  point 
differences  (linear  variation)  in  place  of  the  usual  three  point  difference 
(quadratic  curve  fit) .  This  option  was  applied  to  the  step  region  fi jm  the 
first  station  in  the  step  to  the  second  station  at  the  new  wall  temperature. 

A  data  start  run  for  the  step-down  case  was  made  to  model  the  profiles 
immediately  upstream  of  the  step  exactly.  In  particular,  accurate  total 
temperature  modeling  was  desired. 

(2)  Comments  on  Results 

Profile  and  Cf/2  comparisons  are  presented  in  Figures  66-70  and 

72-76  for  the  step-up  and  r.Lep-down  cases.  Step-down  comparisons  include  both 

zero  start  and  data  start  comparisons.  No  comparisons  are  shown  for  momentum 

* 

thickness  because  the  accuracy  of  these  data  are  being  reviewed  by  TRW.  The 
five  plotted  comparison  profiles  include  the  last  station  before  the  step, 
three  stations  within  the  immediate  region  of  the  step,  and  the  final  station 
which  was  approximately  1  foot  downstream  from  the  step. 

Profile  comparisons  for  both  cases  are  quite  good,  particularly  for  the 
Mach  number  and  velocity.  In  the  totul  temperature  comparison  for  the  step-up 
case  (Figure  67) ,  the  prediction  follows  the  recovery  of  the  profile  well  but 
has  a  larger  "bump"  near  the  wall  in  the  third  and  fourth  stations  shown.  Some 
of  this  is  due  to  the  initial  difference  in  the  profiles  at  the  first  station, 
and  some  may  be  due  to  the  resolution  of  the  data,  that  is,  only  the  open  points 
represent  actual  temperature  probe  data.  The  closest  approach  to  the  wall  was 
0.050  inch;  all  solid  points  shown  inside  this  were  calculated  from  pressure 
probe  data.  A  data  start  prediction  would  have  modeled  the  data  slightly 
better  through  the  step;  however,  the  differences  were  not  regarded  as  serious 
enough  to  warrant  further  computer  solutions. 

The  linear-log  velocity  ratio  profiles  in  Figures  69  and  75  also  reflect 
good  agreement.  Comparisons  for  both  the  first  and  final  profiles  are  shown. 


*  Personal  communication  with  Dr.  Robert  Gran  at  TRW. 
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ction  Coefficient  vs  Streamwise  Location 
p-Up  in  Wall  Temperature  ,  Supersonic 
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Figure  75.  Linear-Log  Velocity  Ratio  Profiles 

TRW,  Step-Down  in  Wall  Temperature,  f jpersonic 
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Figure  76.  Skin  Friction  Coefficient  vs  Streamwise  Location 
TRW,  Step-Down  in  Wall  Temperature,  Supersonic 
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So te  that  stoBtwd  scale*  are  oa  the  ahcisaa.  Clow  rt  app writ  of  the 

3.111  inch  thick  total  ye— art  probe  «at  C.N9  iach,  a  iistaase  well  outside 
the  predicted  laritar  svblajer. 

Fisilly.  tbe  Cj/2  plots  of  Figures  7*  and  S  show  errelTeat  fcweet 
hetwea  the  available  ittia  wall  repoe  Stattse  tabu  aad  predicted  voices.  He 
response  cf  tbe  Cf/2  prediction  to  tbs  step  fol  lew  tbe  voriotioo  io  mil  toe 
per  a  tar*  (Fionss  (S  aed  71>  with  at*  noticeable  la^  iadi  eating  tbe  teak  depeo 
'“sc*  of  tbe  laminar  sab  layer  cm  tbe  oeter  portioo  of  tbe  profile. 

d.  Stanford  Ccepatitess 

fi?  Cj— eats  ob  KW  'spot 

Freest  re  so  aa£  wall  csaditicas  were  enist  aimed  aastaat  withi* 
several  percent  olsasc  tbe  test  section  for  oil  Supscc  ras  used  herein.  Is 
three  of  these  rots  fao  blcaiac,  F  *  I.N2.  sad  step— ag  is  Wcwisj)  .atmgad 
coss'art  strearwise  pressures,  wall  temperatures  and  blowing  aass  flow  rotes 
were  :nrc* .  Fcx  tbe  other  two  ras.  the  octal  reported  voloes  at  each  of  the 
fetr  sfaticos  were  used.  All  predictions  were  started  at  I  *  •  defined  os  the 
Ittdias  edge  of  tbe  test  sectioe.  Tbe  virtual  crisis  resulting  freer  the  use 
of  a  trio  was  located  essentially  at  the  leading  edge  of  tbe  test  sectioc- 

-*2)  Ccoaests  os  Besoits 

Besults  of  Staeted  cteparibxs  are  plotted  is  Figwzes  r»-H. 
Only  velocity  ratio  profile  comparisons  are  iacVed  since  oil  tests  chose*  ere 
essentially  adiabatic  at  the  wall  with  tnpersUre  variations  across  tbe  booad  ■ 
ary  layer  liaitec  to  a  few  percent.  Comparisons  are  presntid  at  each  of  tbe 
soar  stations  reported.  These  ase  at  staticas  he  twees  x  *  1.5  ad  x  =  7-5  feet. 
Station  locations  are-  ac-t  tbe  sate  for  every  test  aad  are  cited  ce  each  plot. 

For  F'qures  85  a*d  8$.  tbe  step  in  blowing  occurs  at  z  *  5.82  fret.  Set  and 
coparisxs  are  pres  rated  in  a  gresp  {Figures  17-?0)  following  tbe  profiles. 

The  four  constant  blowing  runs  are  plotted  together  in  Figures  77  *ad  90. 

Overall  profile  shapes  are  is.  reasonably  good  acreeaeat  at  all  biwiag 
rates  (see  Figures  77-851.  however ,  tbe  oreii'-tioe  indicates  progressively 
higher  values  is  the  aid-range  of  tbe  profile  as  tbe  boundary  layer  proceeds 
dcunstceaa  following  quite  good  agreement  at  the  initial  statioo.  (Doe  to  tiis 
initial  ae.tch i  ng ,  data  start  runs  were  unnecessary.?  This  tread  is  confirmed 
is  the  Se.  plots  of  Figures  87  and  83  which  indicate  slightly  higher  valces 
cf  ?  fez  the  re  as  vi  red  fata.  It  is  el  so  reflected  in  Cf/2  ctwtpa  isnas  is 
Figures  39  and  90.  As  is  well  Juscva,  nine r  errors  in  f  calculation 
can  result  in  eroraous  errors  in  calculated  drag  coefficient  for  flows 
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Figure  78.  Linear-Log  Velocity  Ratio  Profiles 

Stanford,  Simpson  Run  #2367,  Subsonic,  No  Blowing 
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Figure  80.  Linear-Log  Velocity  Ratio  Profiles 

Stanford,  Simpson  Run  #122366,  Subsonic,  F  =  0.002 
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Figure  84.  Linear-Log  Velocity  Ratio  Profiles 

Stanford,  Simpson  Run  #12196C,  Subsonic,  F  =  0.008 
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Figure  86.  Linear-Log  Velocity  Ratio  Profiles 

Stanford,  Simpson  Run  #5867,  Step-Up  in  Blowing 
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with  strong  blowing.  This  holds  true  both  for  the  reduction  of  experimental 
data  and  for  prediction  procedures  such  as  BLIMP.  This  point  will  be  discussed 
in  more  detail  later. 

The  vertical  bands  on  each  data  point  represent  the  reported  uncertainty 
intervals.*  The  differences  between  data  and  prediction  in  te:-ms  of  percent  of 
the  data  value  are  listed  in  Table  XX. 


TABLE  IX 

Cf/2  COMPARISONS  FOR  STANFORD  RUNS 


Simpson  Run 

Difference  in  C^/2 

F  =  0.000 

12%  low 

F  ■  0.002 

25%  low 

F  -=  0.004 

40%  low 

F  =  0.008 

80%  low 

F  =  0.000  *  0.004 

12%  -*■  40%  low 

The  consistency  of  this  pattern  is  completed  by  the  comparison  of  the 
near  wall  region  in  Figures  78,  80,  82,  84,  and  86.  In  every  case,  the  data 
indicate  higher  gradients  at  the  wall  by  roughly  the  percentages  of  Table  VIII. 

The  step  in  blowing  causes  an  immediate  50  percent  decrease  in  predicted 
Cj/2  and  within  the  following  2  feet  (15-20  boundary  layer  thicknesses)  has 
decreased  to  within  10  percent  of  the  value  for  constant  F  =  0.004  blowing  from 
the  leading  edge. 

e .  Jeromin  Comparisons 

(1)  Comments  on  BLIMP  Input 

All  Jeromin  predictions  utilized  the  data  start  option  since  no 
data  were  presented  defining  the  flow  conditions  from  the  nozzle  throat  to  the 
first  measured  profile.  Freestream  conditions  at  each  station  were  based  on 
matching  measured  Mach  numbers;  wall  conditions  (temperature  and  mass  flow) 
were  input  as  reported.  Initial  station  velocity  profiles  were  taken  from 
Squire  (reference  88) .  Temperature  (enthalpy)  profiles  were  computed  from 
equation  (1)  in  Section  II. 3. c. 


*  , 

These  uncertainty  intervals  have  been  described  as  too  narrow  by  Squire, 
reference  87. 
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(2)  Comme  ts  on  Results 

Only  velocity  profile  comparisons  are  presented  because  profile 
data  were  reported  in  that  reduced  form.  For  the  runs  in  question,  total  temp¬ 
erature  profiles  were  not  measured  but  were  calculated  from  equation  (1) . 

This  procedure  was  justified  by  Jeromin  based  on  the  close  agreement  of  pre¬ 
liminary  measurements  with  equation  (1).  Consequently,  since  the  measured 
Mach  number  profiles  and  the  presented  velocity  profiles  are  directly  related; 
no  attempt  was  made  to  convert  back  to  the  Mach  number  form. 

Figures  91-94  and  97-100  contain  the  velocity  profiles  for  the  runs 
2. 5-0.0,  2. 5-1. 2,  3. 6-0.0,  and  3. 6-2.1  in  that  order.  Comparisons  are  pre¬ 
sented  at  four  of  the  five  profile  stations  reported  by  Squire.  These  latter 
profiles  represent  one-half  of  the  number  of  measured  profiles.  Jeromin  reports 
complete  boundary  condition  and  thickness  integral  data  at  all  measured  stations, 
which  covered  a  streamwise  range  from  about  x  *  1.1  to  x  *  1.5  feet.  The  Re0 
plots  in  Figures  95  and  101  include  measured  data  at  all  stations.  A  single 
C^/2  value  obtained  via  the  momentum  integral  method  (by  Jeromin)  is  shown  for 
each  test  in  Figures  96  and  102. 

All  the  velocity  profile  comparisons  indicate  a  recurring  pattern;  the 
BLIMP  profiles  progressively  show  higher  velocity  ratios  in  the  mid-range 
of  the  profile.  By  the  final  profile,  differences  in  velocity  ratios  range 
from  0.05  to  0.07  with  the  maximum  occurring  between  u/u,  =  0.6  to  0.8.  These 
differences  are  only  slightly  greater  for  the  blown  profiles  than  for  the 
unblown . 

The  profiles  reported  by  Squire  are  the  actual  Jeromin  profiles  for 
only  the  Mach  3.6  runs.  Squire  reran  the  Mach  2.5  experiments.  Consequently, 
the  Squire  profiles  and  the  Jeromin  flow  and  wall  conditions  and  computed  pro¬ 
file  parameters  are  related  only  in  a  nominal  sense.  This  appears  to  be  the 
reason  for  the  small  difference  in  initial  data  start  matching  of  the  Re0  values 
for  Mach  2.5  in  the  Figure  95  as  compared  to  the  more  exact  initial  matching  in 
Figure  101  for  Mach  3.6.  The  trends  are  correct,  however,  as  the  Mach  2.5 
profiles  of  Figures  91-94  clearly  indicate  that  the  measured  d0/dx  should  be 
greater  than  that  predicted. 

The  nonlinear  variation  of  Re0  for  Mach  2.5  with  blowing  in  Figure  95 
results  from  the  recorded  but  unintentional  axial  pressure  variation  in  the 
experiment.  The  irregularity  at  the  downstream  end  of  the  test  section  in  the 
Re0  data  in  Figure  101  (Mach  3.6  with  blowing)  is  due  to  apparent  "blowing  off" 
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Fiuure  95.  Momentum  Thickness  Reynolds  Number  vs  Streamwise  Location 
Jeromin  Mach  2.5  Runs 
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Jeromin  Mach  2. 
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Figure  100.  Linear-Log  Volocity  Ratio  Profiles 

Jeromin  Run  3. 6-2.1,  Supersonic,  F  =  0.0021 


JEROMIN  MACH  3.6  RUNS 
NOMINAL  CONDITIONS 


Figure  101.  Momentum  Thickness  Reynolds  Ni"  iiber  vs  Streamwise  Location 
Jeromin  Mach  3 . 6  Runs 
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Figure  102.  Skin  Friction  Coefficients  vs  Streamwise  Location 
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of  the  boundary  layer  as  discussed  by  bcth  Jeromin  and  Squire.  Difficulties 
in  obtaining  pitot  pressure  data  near  the  wall  under  these  conditions  are 
blamed  for  the  unusual  shape  of  the  final  data  profile  (x  ®  1.428  feet)  in 
Figure  99. 

The  Cf/ 2  comparisons  of  Figures  96  and  102  substantiate  the  profile 
and  Re^  comparisons  in  indicating  lower  predicted  values  for  blowing  and 
no-blowing.  The  estimated  uncertainty  intervals  reported  by  Jeromin  are  included. 
The  familiar  date,  start  patterns  are  evident  again.  There  is  an  immediate  change 
to  the  approximate  "equilibrium”  value  desired  by  BLIMP,  followed  by  a  gradual 
variation  caused  by  a  combination  of  axial  pressure  gradients  and,  probably  the 
more  significant  factor,  adjustments  in  outer  profile  shape.  Comparative 
variations  in  the  differences  between  measured  and  predicted  Cf/2  and  d9/dx 
values  relative  to  the  measured  values  are  given  in  Table  X. 

TABLE  X 

VARIATIONS  IN  MOMENTUM  INTEGRAL  TERMS  FOR  JEROMIN  COMPARISONS 


Run 

Cf/2 

d0/dx 

2.5  -  0.0 

18 »  low 

15%  lew 

2.5  -  1.2 

45%  low 

— 

3.5  -  0.0 

15%  low 

28%  low 

3.5  -  2.1 

90%  low 

— 

Values  of  d0/dx  are  included  only  for  the  no-blowing  cases  where,  neglect¬ 
ing  the  effect  of  pressure  gradients,  equation  (38)  applies.  As  discussed  pre¬ 
viously,  the  presence  of  pressure  gradients  does  not  have  a  significant  effect 
on  the  BLIMP  prediction  but  as  noted  by  Jeromin,  it  is  significant  in  the  deter¬ 
mination  of  Cj/2  from  the  momentum  integral  equation.  This  is  particularly  true 
with  blowing  when  both  d8/dx  and  F  are  an  order  of  magnitude  greater  than  Cf/2. 
Considering  the  pressure  gradients  and  the  three-dimensional  effects  reported 
by  Jeromin,  combined  with  the  difficulties  in  accurately  measuring  d0/dx,  it 
appears  likely  that  the  reported  uncertainty  intervals  for  blowing  are  under¬ 
estimated. 

One  additional  prediction  of  C^/2  is  shown  for  Mach  3.6  with  blowing  on 
Figure  102  as  a  solid  curve.  A  zero  start  case  was  run  from  the  nozzle  throat 
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assuming  streamwise  conditions  based  on  descriptions  of  the  nozzle  and  the  phys¬ 
ical  dimensions  of  the  porous  plate  in  references  48  and  88.  The  matching  of  0 
at  the  first  station  was  about  20  percent  off,  but  the  d6/dx  from  that  point  on 
was  essentially  identical  to  that  of  the  data  start  prediction.  The  close  agree¬ 
ment  of  the  solid  and  dashed  curves  again  confirms  the  ability  of  the  data  start 
procedure  to  respond  to  local  wall  conditions. 

Figure  103  shows  a  comparison  of  static  temperatures  at  the  final  station 
for  M  =>  3.6  and  no  blowing.  This  is  presented  to  assess  the  relationships  among 
the  temperature  profile  predicted  by  BLIMP  (dashed  curve) ,  the  temperature  cal¬ 
culated  from  the  Crocco  relation,  using  the  velocity  profile  predicted  by  BLIMP 
(triangles) ,  and  the  measured  Jeromin  profile  (circles) .  The  close  comparison 
between  the  two  BLIMP  determined  temperatures  indicates  that  the  solution  is  in 
close  agreement  with  the  Crocco  relation,  and  thus,  in  agreement  with  Jeromin' s 
observations.  The  difference  noted  between  Jeromin  data  and  BLIMP  predictions 
is  consistent  with  the  difference  in  the  predicted  and  measured  velocity 
profiles  for  x  =  1.428  feet  in  Figure  97. 

2.  RATIONALIZATION  OF  RESULTS 

Very  li*  '  effort  was  expended  to  improve  the  turbulent  model  for  each 
data  set  as  it  was  run,  for  two  reasons.  First,  it  would  be  unwise  to  make 
changes  without  evidence  from  a  number  of  cases  that  a  change  was  called  for. 
Second,  the  contract  for  this  study  did  not  call  for  such  an  optimization. 

Rather,  once  a  model  was  selected,  it  was  to  be  evaluated  for  all  data  sets. 

Now  that  all  the  comparisons  have  been  completed,  however,  it  seems  appropriate 
to  examine  the  results,  suggest  what  improvements  to  make,  and  how  to  make 
them. 

a.  Overview  of  the  Agreement  Between  Experiment  and  Theory 

All  five  of  the  data  sets  which  have  been  used  for  comparisons  here 
include  data  for  unb.lown,  essentially  zero  pressure  gradient  flows.  In  addition, 
the  comparison  with  the  data  of  Wieghardt  and  Tillman  provided  a  baseline  case 
with  which  the  others  can  be  compared.  The  velocity  profile  comparison  was 
given  in  each  case  and  in  general  was  very  good.  Only  the  Lee  case  offered  any 
significant  errors  in  predicted  velocity  profile  shape;  however,  the  validity 
of  the  Lee  data  is  open  to  question.  Mach  number  profiles  for  the  supersonic 
and  hypersonic  data  sets  are  not  predicted  as  well  as  velocity  profiles.  The 
problem  is  generally  one  of  failing  to  predict  an  inflection  point  near  the  mid¬ 
range  of  the  Mach  number  ratio.  In  addition,  M/M  approaches  1.0  with  a  much 
greater  slope  than  the  prediction  shows,  which  of  course  is  related  to  the 
inflection  point  problem. 
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Figure  103.  Static  Temperature  Profiles,  Jeromin,  M  =  3.6,  F  =  0.000 
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Total  temperature  profiles  are  fairly  good  for  MOL  and  TRW  data,  the 
only  cases  where  such  data  were  presented.  The  response  of  the  prediction  to 
the  step  in  wall  temperature  in  the  TRW  data  was  particularly  encouraging.  The 
Cf/2  and  Re^  vs.  streamwise  distance  plots  are  related  through  the  momentum 
integral  equation.  There  seems  to  be  a  general  trend  to  underpredict  C^/2  and 
6  variations  slightly  in  several  of  the  comparisons .  This  trend  is  disappoint¬ 
ing  in  light  of  the  excellent  agreement  with  the  Wieghardt  and  Tillman  data. 

Some  of  this  error  can  be  aliminated  by  optimizing  the  nodal  layout,  as  will  be 
discussed  in  the  next  subsection,  however  it  appears  that  the  prediction  will 
remain  2-3%  low  for  the  better  quality  data,  such  as  Coles  and  Stanford.  Refer¬ 
ring  back  to  the  comparison  between  models  of  Figure  11,  it  :s  clear  that  all 
three  turbulent  models  considered  at  the  beginning  of  this  report  should  give 
about  the  same  results  for  unblown,  low  speed  flows,  with  the  Aerotherm  model 
giving  a  slightly  higher  shear  stress.  Since  the  Cebeci  model  in  particular 
has  been  shown  to  be  very  successful  in  predicting  for  unblown  flows 
(references  68,  77,  and  78),  it  is  hypothesized  that  the  error  observed  here 
is  a  random  one,  ar.d  does  not  indicate  a  trend  associated  with  the  model.  Mere 
will  be  said  of  this  later. 

For  the  blown  flow  data  comparisons,  velocity  profile  predictions  are 
again  quite  good.  In  the  Stanford  data,  there  again  seems  to  be  a  trend  toward 
underpredicting  the  momentum  thickness.  This  results  in  a  gross  underprediction 
of  the  reported  drag  coefficient  for  flows  with  moderate  to  strong  blowing. 

This  result  could  be  anticipated  from  the  on^-dimensional  analysis  comparisons 
of  Figures  24,  25,  and  26.  It  was  clear  that,  for  the  Aerotherm  model, 
both  the  wall  shear  and  the  profile  shape  could  not  be  predicted  simultaneously. 
Thus,  either  the  Simpson  data  are  incorrect,  or  the  Aerotherm  model  should  be 
adjusted  to  fit  it.  This  point  is  discussed  further  in  the  next  subsection. 

b.  Changes  in  the  Turbulent  Model 

For  unblown  flows,  it  has  been  stated  above  that  the  Aerotherm  model 
is  essentially  equivalent  to  other,  apparently  successful  models  and  there  is 
no  reason  why  it  should  not  offer  equivalent  accuracy.  A  significant  improve¬ 
ment  in  accuracy  can  be  made  by  working  with  more  nodes  through  the  boundary 
layer,  as  discussed  in  the  appendix.  As  more  experience  was  gained  in  working 
with  the  code  through  the  course  of  this  contract,  it  became  apparent  that, 
in  addition  to  those  comments  made  in  the  appendix,  the  nodal  distribution  in 
the  transition  region  (figure  104)  rs  of  primary  importance.  Figure  105 
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Figure  104.  Diagram  Showing  Law  of  the  Hall  Nomenclature 

shows  Cf/2  for  a  25-node  model  and  for  two  15-node  runs,  one  with  more  emphasis 
placed  on  the  transition  region  with  correspondingly  less  on  the  wake  region. 
This  weighted  model  is  in  better  agreement  with  the  25-node  run  which  placed 
all  the  extra  nodes  into  the  transition  and  law  of  the  wall  regions.  Figure  105 
also  shows  Cj/2  values  calculated  from  several  popular  theories  at  Re^  *  lit6. 

It  is  apparent  that  there  is  a  certain  amount  of  disagreement  over  the  correct 
drag  coefficient  value  even  fcr  a  simple,  low  speed,  incompressible,  flat 
flat  problem. 

It  is  of  inrerest  to  note  that  the  weighted  15-node  curve  is  nearly  equal 
to  the  standard  15-node  curve  at  the  lower  Reynolds  numbers  but  approaches  the 
25-node  curve  at  the  higher  Reynolds  numbers.  This  is  apparently  due  to  the 
changing  shape  of  the  profile  relative  to  the  fixed  n  distribution.  At  the 
lower  Reynolds  numbers,  the  profiles  are  more  nearly  laminar  with  transition 
occurring  in  the  outer  portion  of  the  n  values.  The  2^-node  model  has  a  suffi¬ 
cient  number  of  nodes  to  model  the  transition  region  at  any  of  the  Reynolds 
numbers  considered,  it  is  obvious  that  the  15-node  model  lacks  this  flexi¬ 
bility,  once  more  emphasizing  the  fact  that  to  obtain  the  most  accurate  pre¬ 
diction,  it  is  necessary  to  evaluate  the  results  with  respect  to  the  chosen 
n  distribution  and  to  select  that  distribution  based  upon  the  streamwise  region 
of  greatest  interest. 

These  comparisons  were  made  after  the  results  described  in  Section  IV 
were  obtained  and  plotted.  Since  the  differences  are  small  percentagewise, 
predictions  were  not  rerun  with  the  improved  n  distributions.  Thus,  all  zero 
blowing  Cf/ 2  predictions  would  appear  to  be  about  5t  low  for  this  reason. 

One  obvious  way  tc  alter  the  turbulent  model  is  to  change  the  numerical 
values  of  the  constants.  As  an  indication  of  how  the  constants  would  affect 
unblown  boundary  layer  predictions,  numerical  experiments  were  performed  for 


153 


BIAS  JUS  PIPE  FLOW 


AFWL-TR- 71-57 


Comparison  of  Predicted  Skin  Friction  Coefficient  for  Zero  Blowin? 
Stanford  Canm  to  Data  and  Fmpirieal  Values 


AFWL-TR-71-57 


the  Coles  M  “  3.7  case.  Alternate  values  of  both  and  y*  were  tried  for  this 
case,  with  the  result  shown  in  Figures  106  and  1C7.  The  small  improvement  in 
Cf/2  prediction  for  both  the  y*  and  k^  modificatior s  was  made,  at  the  particular 
station  selected,  with  a  corresponding  small  improvement  in  profile.  These 
improvements  are,  of  course,  intimately  related  to  the  nodal  distribution 
through  the  transition  region,  therefore  it  is  possible  that  the  results  may 
change  at  much  larger  axial  distances.  It  was  demonstrated  in  Section  III. 3.6 
that  small  changes  in  the  wake  law  eddy  viscosity  constant  also  result  in  small 
but  detectable  prefile  changes  for  unblown  flows.  Eased  on  this  evidence,  it 
appears  that  basic  profile  shapes  and  features  are  unlikely  to  be  changed  with 
modest  changes  in  the  model  constants.  "Fine  tuning'  of  the  turbulent  model 
to  match  drag  or  heat  transfer  data,  for  example ,  may  be  accomplished  with 
small  adjustments  in  these  constants.  Changes  should  be  based  on  more 
comparison  information  than  presented  here,  however.  For  flows  at  low  Reynolds 
number  or  with  high  heat  transfer  rates,  larger  changes  in  these  constants  may 
be  in  order.  The  results  presented  here  will  be  useful  in  estimating  the  results 
of  such  changes. 

For  flows  with  strong  blowing,  the  effects  of  constant  changes  are  altered 
somewhat-  The  wall  law  constants,  y*  and  k^,  have  the  greatest  influence  in 
the  transition  region  of  the  profile,  which  is  much  nearer  the  wall  with  strong 
blowing.  Thus,  there  is  virtually  no  change  in  the  outer  profile  shape  for 
different  wall  law  constants,  as  shown  in  Figure  108.  Significant  differences 
ne.»r  the  wall  do  affect  the  drag  coefficient,  however,  as  seen  in  Figure  109. 
Manipulation  of  the  y  type  constant  is  essentially  the  technique  u,ad  by  both 
Cebeci  and  Bus.inell  and  Beckwith  to  account  for  blowing  in  their  models,  there¬ 
fore  this  path  does  appear  to  be  a  strong  possibility  if  blown  flow  model  changes 
are  indeed  desired.  The  question  of  whether  such  model  changes  are  in  order  is 
addressed  in  the  next  subsection. 

c.  Desirability  of  Turbulent  Model  Changes  for  Blown  Flows 

A  considerable  body  of  data  and  numerous  theories  now  exist  for 
turbulent  boundary  layers  with  injection.  The  discussion,  evaluation,  and  re- 
evaluation,  of  this  data  has  been  a  favorite  topic  in  the  recent  fluid  mechanics 
literature.  Of  particular  interest  for  purposes  of  the  present  discussion  is 
the  drag  coefficient  correction  due  tc  blowing  at  any  given  point  on  a  flat 
plate.  Figure  110,  taken  from  reference  89,  presents  a  number  of 
theoretical  solutions  to  this  problem  for  incompressible  flows,  along  with  a 
few  points  from  Simpson's  data.  The  current  Aerotherm  theory  essentially 
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Unblown  Flat  Plate  Flow 


ure  108.  Effect  of  Wall  Law  Changes  on  Velocity  Profile  in 
Flat  Plate  Flow  With  Blowing 


109.  Effect  of  Wall  Law  Changes  on  Drag  Coefficient  in  Flat 
Plate  Flow  With  Blowing 
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duplicates  the  reference  2  curve.  It  is  clear  that  Simpson's  data  shows  less 
Cf  correction  than  any  of  the  theories  presented.  Figure  111,  taken  from 
reference  87,  illustrates  the  point  that  Simpson's  data  show  less  Cf  correction 
than  that  found  by  other  experimentalists.  Thus,  while  the  Stanford  heat  and 
mass  transfer  apparatus  has  been  very  carefully  constructed  and  operated  by 
competent  researchers,  there  is  not  universal  agreement  that  the  drag  data 
obtained  with  this  apparatus  are  correct.  This  merely  is  a  result  of  the  fact 
that  with  present  measuring  techniques,  the  calculation  of  drag  by  either 
momentum  integral  or  wall  profile  techniques  in  blown  flows  is  not  sufficiently 
precise  to  draw  any  accurate  quantitative  conclusions. 

It  has  been  shown  in  Section  III. 3  that  strong  blowing  is  not  encountered 
in  typical  heatshield  or  nosetip  flight  calculations.  It  is  therefore  concluded 
that  changes  in  the  Aerotherm  wall  law  model  are  not  warranted  until  more 
conclusive  experimental  data  in  the  strong  blowing  region  are  available. 
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SFCTION  V 

CONCLUSIONS  AND  RECOMMENDATIONS 

Results  of  this  research  program  are  summarized  briefly  and  conclusions 
are  drawn  in  Section  V.l.  Recommendations  for  further  work  are  included  in 
Section  V. 2. 

1.  CONCLUSIONS 

The  research  and  development  program  described  in  this  report  has 
uncovered  useful  information  about  the  state  of  the  art  in  turbulent  boundary 
layer  experimentation  and  turbulence  modeling.  It  has  also  answered  many 
questions  about  the  use  of  BLIMP  as  a  prediction  tool.  In  the  literature  survey 
portion  of  this  contract,  it  became  evident  that  no  one  set  of  experimental 
data  is  completely  suitable  for  turbulent  model  studies,  i.e.,  completely  error 
free.  There  is  significant  disagreement  between  various  data  sets  for  even 
very  straightforward  experiments,  such  as  flat  plate  subsonic  flow.  Experimental 
problems  and  potential  errors  are  compounded  for  the  more  difficult  cases,  such 
as  flows  with  blowing. 

Comparisons  of  three  leading  methods  of  modeling  turbulence  in  boundary 
layers  showed  that  the  methods  were  very  similar  in  unblown  flows,  but  contained 
potentially  important  differences  for  flows  with  blowing.  The  one-dimensional 
solution  technique  which  was  used  to  compare  thete  models  appears  to  be 
useful  for  further  development  of  turbulent  models  and  possibly  for  screening 
new  experimental  data. 

The  comparisons  of  predictions  and  data  presented  in  this  report  are 
valuable  in  establishing  the  degree  of  confidence  which  should  be  placed  in 
BLIMP  predictions.  In  general,  it  can  be  concluded  that  the  eddy  viscosity- 
type  model  for  turbulence  is  satisfactory  for  the  type  of  flows  considered  here. 
BLIMP  does  a  good  job  of  predicting  velocity  and  temperature  profiles  for  a 
wide  variety  of  flows  and  boundary  conditions.  Some  improvements  in  profile 
shape  are  possible,  however,  particularly  for  higher  Mach  number  flows.  Of 
particular  interest  in  the  data  profile  shape  is  ar.  inflection  point  in  velocity 
and  Mach  number  profiles  which  occurs  at  supersonic  and  hypersonic  speeds.  Some 
changes  in  the  turbulent  model,  possibly  including  a  variable  turbulent  Prandtl 
number,  may  be  necessary  to  model  this  particular  profile  feature. 
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nr ag  coefficient  predictions  for  many  of  the  no-blowing  cases  considered 
in  this  report  were  slightly  low.  This  was  found  to  be  a  function  of  the 
number  of  nodes  used  and/or  their  spacing  through  the  boundary  layer.  It  can 
be  concluded  that  the  use  of  15  nodes  is  very  near  the  lower  limit  for  accurate 
turbulent  boundary  layer  predictions,  and  that  strong  consideration  should  be 
given  to  up-dimensioning  the  code  to  25  nodes. 

For  boundary  layer  flows  with  blowing,  profiles  were  again  good.  Drag 
coefficients  were  typically  below  the  reported  values  for  both  the  Simpson  and 
the  Jeromin  data.  With  the  Simpson  data,  there  is  considerable  disagreement 
in  the  literature  as  to  whether  the  reported  drag  coefficients  are  correct. 

For  this  reason,  the  rather  poor  agreement  with  the  BLIMP  predictions  could 
only  be  termed  "consistent  but  inconclusive."  Since  strong  blowing  is  not 
typically  encountered  in  heatshield  ablation  problems,  it  is  concluded  that  a 
change  in  the  turbulent  model  to  fit  the  Simpson  or  Jeromin  drag  data  is  not 
justified  without  further  study. 

The  very  large  number  of  computer  runs  necessary  for  the  preparation 
of  accurate  predictions  has  resulted  in  new  information  on  the  use  of  the  BLIMP 
code.  A  technique  for  starting  a  problem  with  a  known  profile  at  the  first 
station  was  developed,  and  the  "data  start"  runs  emphasized  some  interesting 
features  of  the  downstream  solution.  The  general  conclusion  to  be  made  from 
the  data  start  runs  is  that  wall  region  profiles  (and  the  associated  wall  shear) 
approach  the  zero  start  predictions  very  rapidly,  with  the  outer  profile  taking 
somewhat  longer.  This  result  then  provides  some  information  as  to  the  accuracy 
inherent  in  the  usual  technique  of  starting  a  solution  far  upstream  of  the 
region  of  interest,  with  the  expectation  that  starting  profile  errors  will 
die  out  quickly. 

The  comparisons  included  in  this  report  include  demonstrations  of  the 
sensitivity  of  the  predicted  profiles  to  different  numerical  values  of  the 
turbulent  model  constants.  It  is  concluded  that,  for  unblown  boundary  layers, 
small  (< 30%)  changes  in  these  constants  will  not  have  any  important  effects 
on  profile  shapes.  On  the  other  hand,  for  flows  with  blowing,  changes  in  the 
constants  can  have  large  effects  on  the  profiles  vary  near  the  wall,  and 
therefore  affect  drag,  heat  transfer,  etc.  Thus,  changes  in  the  model  constants 
as  a  function  blowing  rate  offeks  a  straightforward  method  of  altering  wall 
parameters  with  this  turbulent  model.  It  may  also  be  concluded  that  the  model 
would  be  sensitive  to  other  changes  such  as  the  substitution  of  tw  for  t  in 
equation  6. 
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2 .  RECOMMENDATIONS 

Many  areas  for  further  investigation  have  become  apparent  during  this 
study.  Perhaps  the  foremost  of  these  is  the  need  for  a  straightforward  con¬ 
tinuation  of  the  kind  of  work  reported  here,  i.e.,  documentation  of  the  validity 
of  the  code  through  comparison  with  experimental  data.  Fine  tuning  of  the 
turbulent  model  for  unblown  flows  should  be  carried  out  in  order  to  make  BLIMP 
the  accurate  and  sophisticated  prediction  tool  which  it  is  intended  to  be. 
Initial  studies  should  concentrate  on  drag  data,  then  be  extended  to  hect  trans¬ 
fer.  Much  of  the  data  screening  and  model  development  work  could  be  carried  out 
most  efficiently  with  a  one-dimensional  code  such  as  the  WALAW  program  described 
in  this  report. 

Once  this  fine  tuning  phase  is  completed,  attention  should  be  given  to 
other  types  of  flows  not  covered  in  the  present  study.  Among  the  many  flow 
regimes  which  merit  attention  are  flows  with  large  heat  transfer,  low  Reynolds 
numbers  (near  transition),  chemically  reacting  flows,  and  flows  in  adverse  pres¬ 
sure  gradient.  All  these  conditions  exist  at  the  surface  of  a  reentry  vehicle, 
where  the  code  is  used  to  predict  the  resulting  boundary  layer.  The  need  for 
verification  is  obvious. 

The  question  of  drag  prediction  in  flows  with  blowing  should  also  be 
resolved,  perhaps  through  comparison  with  wall  heat  transfer  rather  than  drag 
data.  As  a  minimum,  comparisons  with  some  of  the  other  cases  shown  in 
Table  II  should  be  carried  out. 

Finally,  the  need  for  additional  experimental  data  in  all  types  of  flows 
is  apparent.  The  hypersonic  boundary  layer  area  is  of  most  interest  for 
reentry  vehicle  purposes.  Research  programs  aimed  at  the  development  of  new 
instrumentation  for  use  in  blown  or  ablating  boundary  layer  flows  are  particu¬ 
larly  needed  in  order  to  eliminate  the  uncertainties  that  were  brought  out  in 
this  report. 
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APPENDIX 

GENERAL  DISCUSSION  OF  COMPUTER  CODE  SETUP  AND  OPTIONS 


This  appendix  contains  information  on  the  use  of  the  BLIMP  code  for 
the  type  of  problems  encountered  with  the  selected  data  sets.  It  also  presents 
details  on  the  new  entropy  layer  option. 


1.  SELECTION  OF  NODAL  (n)  DISTRIBUTION 

The  BLIMP  solution  procedure  operates  in  the  coordinate  system, 

where  these  quantities  are  defined  as 
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a  (s.n)  grid  system  is  assumed  to  be  superimposed  on  the  boundary  layer  region, 
where  h  is  measured  normal  to  the  wall  and  5  is  measured  parallel  to  it 
(see  Figure  112).  The  boundary  layer  is  divided  into  .,-1  strips  connecting 
N  nodal  points  at  each  £  station.  These  nodal  points  are  designated  by 
where  i  ■  1  *t  the  wall  and  N  at  the  edge  of  the  velocity  boundary  layer.  The 
nodal  system  expands  and  contracts  with  t-.u  boundary  layer  flow,  since  is 
defined  to  be  located  at  the  wall,  and  is  defined  to  be  the  outer  edge  of 
the  boundary  layer.  While  the  (,  numerical  calues  are  calculated  automatically 
by  the  program  from  the  axial  station  dimension  (s) ,  the  numerical  values 
which  *t  uses  are  supplied  directly  as  input.  This  subsection  provides  some 
insight  into  the  selection  of  values  for  proper  program  operation. 

a.  Number  of  Nodes  Requirea 

Since  BLIMP  nolves  a  linear  matrix  of  order  proportional  to  N  (the 
number  of  n  nodes),  the  time  required  to  obtain  a  solution  can  be  expected  to 
be  roughly  proportional  to  N  squared.  Consequently,  it  is  desirable 

*  An  average  of  the  matrix  inversion  (-N*)  and  other  operations  which  are 
proportional  to  N  and  N: . 
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to  determine  the  minimum  number  of  nodes  (and  spacing  of  those  nodes)  which  is 
consistent  with  an  accurate  and  stable  calculation  of  the  boundary  layer.  No 
attempt  was  made  to  evaluate  the  minimum  number  aspect  of  this  question;  however 
on  a  number  of  different  runs  of  Coles,  Stanford  and  others,  two  BLIMP  runs  were 
made  which  were  identical  with  the  exception  of  the  number  of  nodes,  15  being 
used  on  one  and  25  on  ♦•he  other. 

The  results  of  one  of  the  runs  which  are  representative  of  all  those 
made  are  as  follows.  First,  for  runs  of  Coles'  test  #20  with  the  same  8 
axial  stations,  the  number  of  iterations  to  a  solution  at  each  station  was  the 
same  for  either  15  or  25  nodes.  The  time  required  per  iteration  averaged  0.47 
seconds  for  the  15-node  run  compared  to  1.46  seconds  for  the  25-node  run,  a 
ratio  of  0.31.  The  ratio  of  squares  is  0.36;  that  is,  a  25-node  iteration 
took  slightly  longer  than  estimated  by  the  N-squared  proportionality.  The 
differences  in  skin  friction  coefficient,  C^,  and  momentum  thickness,  6,  were 
consistently  about  2.5  percent;  the  15-node  run  having  the  lower  values.  Other 
comparison  runs  indicated  similar  differences  (0  to  5  percent  in  Cf  and  0)  with 
the  15  node  model  giving  consistently  lower  values. 

Figure  113  shows  the  two  velocity  ratio  profiles  on  a  linear-log  scale  with 
the  25-node  run  as  a  solid  line  and  the  15  node  as  circles.  This  evidence  to¬ 
gether  with  that  above  is  judged  sufficient  to  conclude  that  the  15-node  model 
represents  a  worthwhile  saving  in  computer  time  while  maintaining  satisfactory 
accuracy.  Consequently,  this  model  was  used  extensively  nr  making  the  final 
BLIMP  predictions.  Due  to  the  smaller  number  of  nodes  available,  however,  care¬ 
ful  judgement  had  to  be  exercised  in  the  choice  of  the  nodal  distribution. 

b.  Distribution  of  Nodes 

Figure  114  depicts  a  typical  turbulent  velocity  profile  and 
also  che  variation  of  the  first  derivative  of  velocity  through  the  boundary 
layer.  The  velocity  gradient  typically  decreases  three  to  four  orders  of  magni¬ 
tude  between  the  wall  value  and  the  value  at  u/ue  ■  0.9.  Thus,  th?  a  priori 
selection  of  a  nodal  spacing  to  "curvefit"  these  variations  with  ten  or  fifteen 
discrete  values  is  a  difficult  problem. 
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Figure  114.  Typical  Velocity  and  Velocity  Gradient  Profiles 
in  a  Turbulent  Boundary  Layer 

The  BLIMP  manual  (reference  90)  suggests  that  the  nodal  spacing  should 
be  such  that  each  successive  value  of  n  not  exceed  the  previous  value  by  much 
more  than  a  factor  of  2.  This  guideline  has  proven  to  be  generally  v?  lid  but 
should  be  evaluated  relative  to  each  particular  type  of  profile.  As  profiles 
become  distorted,  as  in  the  case  of  blowing,  it  may  be  necessary  to  warp  the 
distribution  as  well  to  be  certain  that  regions  of  high  gradients  are  repre¬ 
sented  adequately.  For  curvefitting  purposes,  it  is  undesirable  to  have  a 
change  of  greater  than  0.1  in  the  velocity  ratio  between  any  adjacent  nodes. 

If  experimental  data  are  available,  a  quick  study  of  reported  velocity  profiles 
will  enable  the  user  to  select  a  satisfactory  distribution.  In  the  event  data 
are  unavailable,  a  short  computer  run  limited  to  a  few  stations  can  confirm 
the  adequacy  of  the  selected  distribution  or  indicate  necessary  changes. 

Perhaps  the  most  important  part  of  the  distribution  is  that  nearest  the 
wall.  Since  BLIMP  calculates  Cf  from  the  wall  velocity  gradient,  and  since 
the  gradient  at  the  wall  is  taken  as  the  first  derivative  of  the  first  spline 
fit  quadratic  evaluated  at  y  »  0  (see  reference  1  ) ,  it  is  essential  that  the 
first  several  nodal  points  be  located  within  the  laminar  sublayer.  As  a  rule  of 
thumb,  at  least  the  first  two  points  away  from  the  wall  should  have  velocity 
ratios  less  than  0.1.  This  together  with  the  general  spacing  guideline  above 
should  result  in  satisfactory  prediction  of  the  various  wall  parameters. 

Another  less  obvious  region  of  concern  is  the  outer  edge  of  the  boundary 
layer.  The  properties  of  the  quadratic  and  cubic  spline-fitu  (the  outer  two 
points  are  join*  i  by  a  cubic)  are  such  that  injudicious  nodal  spacing  near  the 
outer  edge  can  cause  an  oscillatory  overshoot  of  the  edge  value  of  1.000  as 
shown  in  Figure  IIS.  In  this  example  the  second  to  last  node  has  been 
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Figure  115.  Overshoot  of  the  Spline  Fit  Procedure 


chosen  as  the  fixed  node  (generally  defined  as  u/u  *  0.95  and  n  3  1.000).  The 

e 

next  to  last  node  may  have  been  chosen  too  close  to  the  fixed  node  and/or  the 
last  node  too  far  from  the  fixed  node.  Since  the  cubic  is  constrained  to 
u/u  =  1.00  and  3(u/u  )/3n  ■  0.00  at  the  final  point  (in  the  absence  of  an 
entropy  gradient),  it  is  possible  to  generate  an  overshoot.  Such  an  occurrence 
is  especially  troublesome  since  integral  properties  are  computed  by  an  exact 
integration  of  the  quadratic  and  cubic  curve  segments.  If  the  differences  in 
n"  are  large,  e\en  small  excursions  in  u/ufi  may  lead  to  large  errors  in  the 
various  thickness  integrals.  Avoidance  of  this  problem  is  again  possible 
by  observing  simple  guidelines. 

First,  if  a  data  profile  is  available,  the  spacing  of  the  final  nodes 
can  be  approximated  well  enough.  For  subsonic  flc>% .  *7  and  y  are  directly 
proportional.  Supersonic  flow  requires  somewhat  smaller  n  spacing  due  to  the 
(generally)  decreasing  density.  Flows  with  blowing  tend  to  approach  the  edge 
condition  with  higher  gradients  (du/dy)  and  thus  require  smaller  spacing  com¬ 
pared  to  unblown  flows.  If  a  profile  is  not  available,  a  trial  run  may  be 
necessary  using  the  generalised  appro..  v  described  below.  Special  attention 
should  be  given  to  the  values  of  FP  (F  PRIME  ■  u/u  )  and  to  FPP  (F  DOUBLE 

C 

PRIME  ■  3 (u/u  ) /3ni  at  the  last  few  points,  particularly  those  points  between 
the  fixed  point  and  the  final  point.  If  an,  of  the  points  has  a  velocity 
ratio  very  near  or  exceeding  1.00  and/or  if  there  is  an  inflection  or  change 
of  sign  in  the  derivative,  an  overshoot  may  have  occurred.  It  is  always  help¬ 
ful  to  plot  these  two  values  over  the  final  few  points  on  a  linear  seal*'  to  be 
sure. 
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During  the  course  of  this  investigation  the  following  choices  for  the 
~  distribution  were  made  and  were  found  to  provide  satisfactory  results.  The 
selection  of  which  node  should  be  the  fixed  ~  node,  the  “  value  for  that  node 
and  the  velocity  ratio  at  that  node  is  arbitrary.  Choosing  the  value  of  o'  as 
1.00  makes  for  ease  in  ratioing  to  other  values  and  is  traditionally  accepted, 
as  is  the  choice  of  0.95  for  the  velocity  ratio  (turbulent  flow  only).  Both 
have  been  used  exclusively  for  final  BLIMP  runs  presented  in  this  report.  The 
number  of  the  fixed  node  has  been  chosen  as  the  13th  out  of  15  total  nodes. 

In  some  earlier  runs,  the  12th  was  fixed;  however,  this  increases  the  possi¬ 
bility  and  the  amplitude  of  overshoot  due  to  a  poor  choice  of  spacing  and 
additionally,  places  more  nodes  than  necessary  in  a  region  of  the  profile 
that  does  not  require  as  much  detail.  The  next  to  last  node  should  be  chosen 
so  as  to  result  in  a  velocity  ratio  of  0.980  +  0.005.  This  will  minimize  the 
possibility  of  overshoot.  Once  again,  an  available  profile  is  the  best  guide. 
The  values  in  Table  XI  which  ere  used  for  this  contract  may  be  used  as  a  start¬ 
ing  point  if  profiles  are  not  available. 

TABLE  XI 

TYPICAL  n  DISTRIBUTIONS  IN  THE  OUTER  WAKE  REGION 


Type  of  Flow 

' 

14  th  n'  Point 

15th  n  Point 

Range  of 
Values 

Mean 

Value 

Range  of 
Values 

Mean 

Value 

Subsonic,  no  blowing 

1.3 

1.3 

1.8 

1.8 

Subsonic,  blowing 

1.2  +1.3 

1.25 

1.5  +  1.7 

1.6 

Supersonic,  no  blowing 

1.3  +  1.5 

1.4 

1.8  +  2.7 

2.1 

Supersonic,  blowing 

1.15  +  1.2 

1.2 

1.67  +  1.75 

1.7 

Hypersonic 

1.5  +  1.7 

1.6 

2.5  +  3.0 

2.'i 
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One  final  point  should  be  considered.  As  the  solution  procedure  pro¬ 
gresses  axially  in  the  flow  direction,  the  rT  grid  system  is  "stretched"  in 
order  to  remain  fixed  to  the  outer  edge  of  the  growing  boundary  layer.  Since 
the  laminar  sublayer  does  not  grow  as  rapidly  as  the  turbulent  outer  flow,  the 
grid  points  nearest  the  wall  may  eventually  be  stretched  out  of  tha  laminar 
region  altogether,  thereby  giving  inaccurate  wall  gradient  information.  It  is 
necessary  therefore  to  check  the  results  at  all  important  stations  to  be  certain 
that  the  n"  distribution  is  sufficient,  especially  at  the  wall.  Restarting  the 
program  at  some  intermediate  body  station  is  a  possibility  for  very  long  running 
lengths. 

2.  THE  "DATA  START"  PROCEDURE 

In  the  classical  boundary  layer  problem,  the  initial  and  boundary 
conditions  ar<‘  such  that  a  known  profile  (or  profiles)  is  provided  at  an  up¬ 
stream  station,  and  sufficient  edge  and  wall  conditions  are  provided  along  the 
flow  direction.  Profiles  are  then  found  at  body  stations  of  interest  based  on 
the  boundary  conditions  and  upstream  profiles.  In  typical  AFWL  applications  of 
the  BLIMP  code,  however,  upstream  profile  information  is  not  available,  and  the 
program  has  been  written  to  start  the  solution  based  on  a  similar  solution  pro¬ 
file  at  the  first  station  (usually  close  to  X  o  0).  When  started  in  this 
manner  sufficiently  far  upstream  cf  the  region  of  interest,  the  errors  resulting 
from  the  approximate  profile  die  out  and  accurate  solution”  are  obtained  at 
downstream  locations.  For  most  of  the  data  sets  chosen,  edge  and  wall  infor¬ 
mation  were  available  in  sufficient  detail  and  accuracy  to  enable  starting  BLIMP 
in  the  conventional  manner,  that  is,  with  a  similar  solution  profile  at  an  axial 
station  slightly  greater  than  zero  (chosen  here  as  0.01  feet).  In  the  super¬ 
sonic  and  hypersonic  tunnel  flows  where  the  test  surface  forms  one  side  of  the 
nozzle,  conditions  through  the  acceleration  section  are  often  not  defined.  Such 
was  the  case  with  Jeromin,  Lee,  et  al.,  and  to  some  degree  with  Brott,  et  al. 

An  alternate  procedure  to  enable  comparisons  to  be  made  is  to  start  DLIMP  in  the 
more  classical  manner,  i.e.,  at  the  first  reported  data  profile  station  with  the 
measured  velocity  and  temperature  (enthalpy)  profiles.  The  correct  imple¬ 
mentation  of  this  starting  procedure  requires  some  detailed  consideration, 
however 

The  BLIMP  code  instructions  as  described  in  reference  90  do  indicate  that 
the  program  can  be  started  with  an  "input"  profile  as  the  solution  at  the  first 
station.  This  starting  option  was  intended  primarily  for  restarting  purposes, 
wherein  an  actual  BLIMP  solution  is  used  as  input  to  the  code  for  further  calcu¬ 
lations.  As  an  actual  solution,  tni*  input  profile  satisfies  the  spline  fit 
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quadratic  and  cubic  relations  between  nodes  which  are  inherent  in  the  integral 
matrix  formulation.  Thus,  to  use  this  starting  option  successfully  with  actual 
experimental  profiles  at  the  first  station,  the  spline  fit  requirement  must 
also  be  satisfied.  For  the  "data  start"  cases  presented  in  this  report,  another 
computer  program*  was  used  which  accepted  the  measured  velocity  and  temperature 
profiles  in  physical  coordinates,  converted  them  to  the  n  coordinate  system, 
performed  a  least  squares  spline-type  curve  fit  for  15  preselected  “  nodal  values, 
and  printed  out  the  stands' d  BLIMP  "restart"  information.  This  restart  infor¬ 
mation  was  then  provided  to  BLIMP  as  the  first  station  profile.  All  data  start 
runs  reported  herein  include  this  technique  and  in  all  cases  convergence  was 
achieved  at  the  second  station  ir.  a  normal  manner,  that  is,  in  from  3  to  8 
iterations  with  an  average  of  only  four  iterations  . 

A  few  comments  regarding  the  preparation  of  data  and  the  resulting  per¬ 
formance  of  the  data  start  procedure  arc  in  order  at  this  point  to  establish 
the  degree  of  accuracy  which  has  beer:  attained  in  the  use  of  this  technique. 

Since  none  of  the  profiles  of  presen  interest  are  sufficiently  defined  near  the 
wall,  it  is  necessary  to  generate  additional  points  for  input  to  the  least 
squares  curve  fit  program.  Figure  116  shows  a  typical  case:  the  data  of  Jeromin 
run  2. 5-1. 2  (Mach  2.5,  F  *  0.0012).  Only  the  portion  of  the  profile  near  the 
wall  is  shown;  the  circles  represent  the  reported  data  and  the  dashed  line  rep¬ 
resents  the  value  of  d(u/ue)/dy  at  the  wall  based  on  the  reported  value  of  Cf. 

The  solid  curve  represents  an  "eyeball "curve  fit  from  which  additional  points 
were  chosen.  The  necessity  of  this  is  apparent  considering  that  at  least  two 
values  of  7  not  including  n  *  0.0  should  be  chosen  with  u/ug  <  0.1  and  about 
3  more  chosen  in  the  remaining  interval  up  to  the  first  reported  data  point  away 
from  the  wall.  Any  other  regions  which  may  not  contain  a  sufficient  number  of 
data  points  may  be  "filled  in"  in  a  similar  manner.  The  resulting  supplemented 
experimental  data  constitute  the  required  input  profile.  The  temperature  pro¬ 
file  is  obtained  in  the  same  manner  except  in  Jerowin’s  cases,  wherein  tempera¬ 
ture  was  defined  by  the  velocity  profile. 

It  is  also  of  interest  to  examine  the  solution  r.t  nearby  stations  once 
BLIMP  accepts  the  data  start  profile  and  proceeds  on  downstream.  Figure  117 
shows  the  progression  in  profile  shape  for  the  Jeromin  case  mentioned  above. 

The  solid  curve  with  circles  indicates  the  input  velocity  profile.  The  subse¬ 
quent  change  in  the  solution  for  each  nodal  point  is  shown  by  the  various 
symbols  noted  in  the  legend  (including  the  station  location).  These  solutions 

*  This  program  was  on  hand  at  Aerotherm  and  was  not  developed  under  this 
contract. 
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Figure  116.  Definition  of  Near  Wall  Velocity  Profile  for 
Data  Start  Input 
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are  also  connected  by  lines  indicating  the  loci  of  the  nodal  point  solutions. 

The  nodal  points  of  the  final  calculated  profile  are  connected  by  the  dashed 
curve.  From  comparisons  made  on  other  cases,  this  dashed  curve  would  be 
quite  similar  to  the  profile  shape  BLIMP  would  predict  based  on  a  zero  start. 

At  the  second  axial  station  (triangle  symbols) ,  there  is  little  if  any 
change  over  the  bulk  of  the  profile;  however,  near  the  wall  major  adjustments 
have  occurred.  This  is  due  to  a  difference  between  the  input  (d(u/ue)/dy 
at  the  wall)  and  the  value  BLIMP  would  predict  given  the  same  free  stream  and 
wall  conditions.  This  adjustment  takes  place  over  very  small  axial  distances. 

It  appears  t-,  be  identical  to  adjustments  which  are  made  in  response  to  steps 
in  wall  conditions  such  as  temperature  and  blowing  rate.  Although  the  magnitude 
of  the  adjustment  is  large,  it  is  limited  to  a  very  thin  layer  near  .he  wall.  It 
therefore  has  a  negligible  effect  upon  thickness  integrals  which  change  in  a 
smooth  and  continuous  manner  (see  Section  IV) .  As  the  solution  continues  on 
downstream,  an  adjustment  is  made  to  the  outer  portion  of  the  law  of  the  wall  re¬ 
gion  while  the  overall  "normal"  growth  of  the  boundary  layer  is  reflected  by  a 
steady  thickening  out  in  the  wake  region.  By  the  final  dashed  profile  shown,  the 
adjustments  to  the  shape  appear  to  be  complete,  and  ordinary  boundary  layer  growth 
accounts  for  any  shifting  of  the  points.  A  number  of  examples  of  the  manifesta¬ 
tions  of  this  behavior  are  apparent  in  the  graphical  presentations  in  Section  IV. 
In  general,  "recovery"  from  the  data  start  profile  is  (1)  essentially  immediate 
for  wall  properties  such  as  Cf,  (2)  for  all  practical  purposes,  unaffected  by 
the  starting  values  of  the  various  integral  parameters  and  (3)  slow  for  the 
overall  profile,  taking  on  the  order  of  10  boundary  layer  thicknesses.  This 
latter  "recovering"  simply  implies  that,  given  the  existing  conditions  at  the 
first  reported  profile  station,  BLIMP  would  not  ha*,  predicted  that  profile, 
the  difference  in  shape  being  similar  to  that  between  the  solid  and  dashed  pro¬ 
files  in  Figure  117. 

3.  HOMOGENEOUS  FLOW  CONSIDERATIONS 

All  of  the  cases  considered  in  this  turbulent  model  study  involved  only 
air  as  a  working  fluid,  both  as  the  main  stream  gas  and  as  the  transpirant. 

Also,  the  temperature  range  under  consideration  was  low  enough  that  no  chemical 
reactions  would  take  place.  Thus,  it  is  possible  to  realize  some  economies  in 
the  operation  of  the  BLIMP  program  in  this  limited  thermochemical  regime.  The 
program  was  modified  to  operate  in  a  homogeneous  flow  mode  by  accepting  a 
single  species  in  the  thermochemical  data  deck.  This  species  is  treated  as  an 
element  (named  "cold  air");  it  is  given  its  own  fictitious  atomic  number  and 
the  usual  set  of  thermochemiral  curvefit  constants.  Accurate  transport 
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properties  can  also  be  realized  in  this  mode  of  operation  by  'eading  in  correct 
diffusion  factor  data  under  Group  12  of  the  program  input. 

4.  ENTROPY  LAYER  OPTION 

A  new  feature  of  the  BLIMP  code  is  the  entropy  layer  option,  which  offers 
a  direct  coupling  of  entropy  gradients  in  the  inviscid  flow  with  the  boundary 
layer  edge  conditions.  Inviscid  flow  entropy  can  be  determined  directly  as  a 
function  of  the  dimensional  stream  function.  Further  this  functional 
relation  can,  for  the  las »  of  a  typical  reentry  vehicle,  be  determined  from  the 
shock  shape.  The  entropy  is  established  from  the  shock  angle,  and  the  stream 
function  from  the  simple  integration  of  a  uniform  freestream  flow. . 

If  the  boundary  layer  'swallows"  a  sufficient  flow  mass  to  result  in 
sizeable  entropy  variations  within  the  swallowed  mass,  it  is  necessary  to 
account  for  this  phenomena  adequately.  An  iterative  mass  balancing  procedure 
is  used  to  establish  the  correct  edge  entropy.  Given  an  initial  estimate  of 
the  edge  condition,  a  solution  is  generated.  The  edge  stream  function  from 
this  solution  demands  a  certain  edge  entropy,  which  is  then  compared  with  the 
edge  entropy  calculated  from  the  edge  pressure  and  enthalpy.  These  two  entro¬ 
pies  may  be  different,  in  which  case  the  estimate  of  edge  conditions  is  altered 
and  a  new  solution  is  obtained.  This  iteration  procedure  is  an  integral  part 
of  the  existing  solution  iteration  procedure,  and  therefore  does  not  generally 
add  to  the  number  of  iterations  required  or  the  program  run  time. 

The  fact  that  the  BLIMP  code  is  written  in  terms  of  normalized  dependent 
variables  adds  some  confusion  to  the  entropy  layer  solution.  It  is  well  known 
that  the  existence  of  an  entropy  gradient  at  the  edge  of  the  boundary  layer  also 
requires  a  velocity  gradient,  i.e., 


3s 

W 


(42) 


Since  velocity  varies  in  the  inviscid  flow  at  the  edge  of  the  boundary  la/er, 
the  definition  of  a  u^  value  to  be  used  in  forming  the  diw  -isionless  velocity 
ratio  u/u#  is  difficult.  The  problem  was  solved  by  using  a  reference  velocity, 
ur,  defined  by  an  isentroplc  expansion  from  stagnation  conditions.  The  formula¬ 
tion  of  the  equations  as  carried  out  in  reference  1  remains  valid  with  the 
new  stipulation  thet 
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fj_  =  u_ 

aH  *  ur  (43) 


The  pressure  gradient  parameter  $  enters  the  nondimens iona)  formulation  of  the 
momentum  equation  in  the  same  way,  since  the  Bernoulli  equation  holds  for  the 
isentropic  expansion 


1  3P 


p  Jx 


9ur 
ur  Tx~ 


(44) 


and  we  will  define 


3  in  u 

e  H  2  nr? 


(45) 


The  pressure  gradient  term  in  the  momentum  equation  remains  essentially 
unaffected  as 


In  the  BLIMP  solution  procedure,  the  value  of  ug/ur  =  f^/a H 
is  a  variable  at  each  station  for  entropy  layer  flows.  It  is  determined  from 
straightforward  energy  relationships  for  the  reference  and  actual  inviscid 
expansions  around  the  body.  As  with  the  isentropic  edge  condition,  the 
selection  of  the  solution  domain  (i.e.,  the  selection  of  the  maximum  value 
of  FT)  is  arbitrary.  Conventional  techniques  for  interpreting  the  results  in 
this  solution  domain  must  be  re-examined,  however.  For  example,  at  any  given 
body  station,  two  individual  computer  runs  with  different  choices  for  the 
numerical  value  of  n"N  would  result  in  two  different  edge  velocity  values. 
Different  values  cf  u  /u  would  also  occur.  Both  solutions  are  correct, 

©  Jt 

however,  since  a  velocity  gradient  should  exist  in  an  entropy  layer  region. 

This  velocity  gradient  will  project  one  edge  state  to  the  other  and  the  choice 
of  the  nN  value  merely  determines  how  far  into  the  inviscid  region  the  boundary 
layer  solution  will  extnd. 

One  other  interesting  feature  of  the  entropy  layer  operation  of  the 
program  is  the  definition  of  the  coordinate  stretching  parameter,  a(J.  Formerly, 
coordinate  stretching  was  accomplished  by  constraining  some  arbitrary  point 
near  the  boundary  layer  edge,  ?c#  to  have  a  specified  velocity  ratio,  C,  neat 
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(but  something  less  than)  the  edge  value  (reference  1) .  With  an  edge  velocity 
gradiertt,  this  constraint  has  been  modified  to  deal  with  the  straight  line 
extrapolation  of  the  edge  velocity,  rather  than  the  ratio  itself,  as  shown 
in  Figure  118.  The  selection  of  kappa  and  CBAP.  input  values  (Group  4,  Card  3 
reference  90)  is  not  changed. 
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Figure  118.  a  Constraint  With  Entropy  Layer  Flow 
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